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EDITORIAI 


SALUTATIONS TO MAURICE JOLY, HAROLD WAYLAND AND THE NEW EDITORS 





It is with great regret that Maur : y and Harold iyland, due to 


their retirement from the institutions they sre associated with for many 


years, found it necessary to terminate their appointments as Editors of 


RHEOLOGY. The Publishers and I should like to express our gratitude t« 


of them for their many contributions to our Journal. 


Maurice Joly has been an Editor of BIORHEOLOGY since the beginning 
journal in 1962. He also acted as Acting Co-Editor-in-Chief from June 
September 1975 when he kindly took over the duties of George W. Scott Blairz 
who was incapacitated due F accident ind eye operations During this t 
he has been most helpful in making important editorial deci ms and 
editing with me the Proceedings of the Second International 
rheology, held at the Weizmann Institute of Science in Reho 
congress contributions were published in several issues of o journal. 
Editor, he has been most helpful during the nineteen years of his associ 
with our journal, including the past three years when BIORHEOLOGY became 
journal for rapid publication, ; > 7 when BIORHEOLOGY b 


a bi-monthly publication. 


The special issue for the Third International Congress of Biorheology 
Volume 15, Numbers 5/6, 1978 was dedicated to Maurice Joly. [t contains 
port I wrote on the Poiseuille Award dedicated to him by the International 
ciety of Biorheology, Syoten Oka's presentation address at the Poiseuille 
Medal Award Ceremony and Maurice yly's Poiseuille Award Lecture. The double 
issue also listed nearly thirty main publications by Maurice Joly and associ- 
ates. Maurice Joly's counsel in many editorial decisions, in particular dur- 
ing the past few years, has been most helpful to me and I should like to take 
this opportunity to thank him for all he did for our journal 


city as member of the Council of our Society. 


Harold Wayland served well as an Editor the past three 


regret that he could not complete his term of 


The Readers of BIORHEOLOGY will join me very best wishes t 
Joly and Harold Wayland for good health |; é for any scientifix« 
butions they may make during theirz etirement. Our Readers will be interested 
to learn that both intend to continue, as much as possible, with their scien- 


tific research in biorheology and other are: of interest to them. 


[It is a tribute to Maurice Joly, Harold Wayland and our other Editors who 
have served or are serving that BIORHEOLOGY, the official journal of the In- 
ternational Society of Biorheology, is the leading journal in all fields of 
our growing science of biorheology and that our journal enjoys a high reputa 


tion in the scientific community. 


It is with gre: pleasure that the Publishers and I welcome our new Edi- 
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tors, J.C. Healy of Saint-Etienne, France and Mitchell Litt of Philadelphia, 


Pennsylvania, U.S.A., who have taken on the duties of Maurice Joly and Harold 
Wayland. 

[I am confident that our Editors and Authors will continue to maintain the 
high standing of BIORHEOLOGY as a significant international scientific jour- 
nal. 

The Publishers and I extend our best wishes for 1980 to the Readers of 
BIORHEOLOGY. 

Alfred L. Copley 


Editor-in-Chief 
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EDITORIAL 


ON SUBMISSION OF COMMUNICATIONS FOR PUBLICATION 





The Publishers and I must again ask the Authors to foll 
Instructions to Authors, printed on the inside back cover of ach issue 
BIORHEOLOGY. There have been delays of publication, because a communication 
accepted for publication on scientific merits had to be returned to the Au- 
thors for corrections and retyping, because of errors esentation. Actu- 
ally, such a technically faulty communication should never have been conside- 
red for review on scientific merits by the Editor, to whom it was submitted. 
Although we have cautioned our Editors, some of them or their 
sistants have overlooked the requirements, as given in the Instructions 


thors. 


Some errors in a manuscript, submitted for camera-ready printing, occur 
frequently and, therefore, I should point out to them in order to facilitate 
publication as rapidly as possible, after a Paper or Brief Communication has 


been accepted by an Editor. They concern the following points: 


l. Preparation or laysheets should be used as much as possible. In 
case these sheets cannot be quickly secured by Authors, eager to 
get their work published without delay, the limitation in size of 
the typing area should be followed as described in the Instruc- 
tions to Authors. The typing area on each sheet should be uti- 
lized as fully as possible. 


Key Words are always necessary and must be placed as a footnote 
after the wording "Key Words:" on the Title Page and not on a se- 
parate sheet. 





Running Titles of not more than 45 typewriter spaces must be pro- 
vided by the Author on a separate sheet or written in blue pencil 
on the top of the Title Page. 


Figures 


Space of Figures on the preparation or laysheets must be pro- 
vided, but not more than is necessary, and the legends will 
have to be placed next to or below the Figures. 


The Figures should be clear either on a glossy photograph or 
in the original. They can never be smaller than the requi- 
red size, but can be larger. 


The lettering for each Figure has to be li > enough to be 
easily readible after the 20% reduction whi each manuscript 
sheet is subjected to. 


On the back of each Figure, the names of the Authors and the 
number of the Figure should be placed lightly in pencil. 


Figures cannot be placed in the text under Discussion. 


The References must be followed in accordance to our Instructions 
to Authors. Surnames must be capitalized, the name of the jour- 
nal and its volume must be underlined and in the case of a book, 
its title has to be underlined. If the book has particular Edi- 


tors, their names must be provided in accordance with the In- 


structions. The last page number of a reference must be given 
together with the first. 
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ing is no longer permissible. Each 
need t« i le-spaced. Each paragraph should be- 
in indent : nd should be separated from the fol- 
raph by ible pacing. In this way a paragraph 


letails are given in the Instructions 
ommunication exceeds the permissi- 

h are free of charge, the Authors are 

itement signed by them or their insti- 
é harge of dollars per every 
he Publishers. 


forms o the Authors, when a commu- 
publication. One of these forms, 
Order Form, notifies the Author(s) 
ommun i ti ’ been accepted, while the other is 
ight Form. Author must complete both forms and 
with l t the Executive Editorial Office (EEO) 
the Publi e1 The EEO is under my direction and 
t the at of Biorheology, Polytechnic In- 
Street, Brooklyn, NY 11201, USA. 
ilso be completed, whether or not 


requested. 
-d fo corrections of errors which are usu- 
there are other points which are given un- 
ich must also be strictly adhere to by the 
issistant must, therefore, ac- 
Instructions » Auth before the typing of 
After the tific : of acceptance of 
ithors must act also in accordance with the 


emphasized in the Editorial. 


only the last revision of June ‘ our 


in each issue beginning with Volume 16, Num- 


Alfred L. Copley 


Editor-in-Chief 
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a blood transfusion center anywhere must have a reasearch laboratory of hemo- 
rheology, as well as a hemorheology laboratory for routine testing of blood 
samples. In any hospital in the world, where blood transfusions are performed, 
the blood is tested for blood groups and other factors. It is very clear to 
everyone, including the lay public, that these are immunological tests for 
compatibility. However, it did not occur thus far to the medical profession 
in general and especially to many blood transfusionists that what they are 
dealing with are hemorheological phenomena on an immunological basis, namely, 
agglutination or the irreversible clumping of red blood cells, and hemolysis. 
It is not yet general knowledge, what we all here assembled take for granted, 
that there can be hemorheological phenomena on a non-immunological basis, 
which can lead following blood transfusion to disturbances in the circulation. 
As I pointed out to my friend van Loghem, there is a need for hemorheological 
studies and routine testing in a blood transfusion center, such as the one in 
Holland, and I reasserted my belief in such a necessity, when I was introduced 


to the new director of this large institution in Amsterdam. 


I know that such hemorheology laboratories do not exist, for instance, as 
part of the New York Blood Center in the city, where I am residing, nor with 
the Blood Service of the American Red Cross in Washington, D.C. To my know- 
ledge, the only blood transfusion center in the world, where hemorheological 
testing and research in hemorheology are being made, sre at the Centre Re- 
gional de Transfusion Sanguine et d'Hematologie. I am glad to say that both 
fundamental and clinical hemorheological research is being pursued here vigo- 
rously and with great enthusiasm. I, therefore, salute Professors Streiff, 
Larcan and Stoltz who had the foresight to include the practice of hemorheolo- 
gy as an integral activity of the transfusion center here at Nancy. The ques- 
tion may be asked, what kind of routine hemorheological testing can be done in 
a blood transfusion center? As I have been told yesterday by Dr. Stoltz, the 
following routine procedures are being employed here at the transfusion cen- 
ter: measurements of viscosity and viscoelasticity of blood and plasma, dif- 
ferent filtration tests, and also electrophoresis tests red blood cells and 
platelets. Since many hemorheological esearch studies >» being made here, 
it can be expected tha as a result, several tests employed in these studies 
will then be applied as routine hemorheoclogical tests for blood transfusion. 
It goes without saying that for routine hemorheological testing, measuring 
techniques and methods developed in othe laboratories should be included as 
far as i feasible or practical, but deems necessary in preventive medicine 
and for clinical considerations. 

should like to emphasize that the ar >» of hemorheology not mere- 


ly imperative in any blood transfusion center, but also in any hospital and 


clinical laboratory, as part of the armamentarium in providing needed informa- 


tion to physicians and surgeons everywhere, for the betterment of medical 


care." 


Alfred L. Copley 
Laboratory of Biorheology 
Polytechnic Institute of New York 


Brooklyn, New York 11201, U.S.A. 
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IN VIVO AND IN VITRO, - A COMPARISON BETWEEN DIFFERENT 
PLASMA SUBSTITUTES 
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Department of Surgery I, University of Gothenburg, Sahlgren’s 
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ABSTRACT The viscosity of whole blood (WB), plasma (P), and packed cells (P¢ 
measured in a Wells- Brookfield cone plate viscometer on blood from dogs 
undergoing surgical shock and hemodilution with different plasma substi- 
tutes and on human ACD-blood diluted in vitro to comparable concentrations 
of colloids, Erythrocyte sedimentation rate (ESR) and microscopic red cell 
aggregate count were used as in vitro indices of red cell aggregation, Du- 
ring three hours of intestinal shock WB viscosity increased by 20 %, hema 
tocrit increased from 50.9 to 59,1 %, P viscosity was unchanged and P< 
viscosity increased by 6 %, When compared at equal hematocrit drop, Ring- 
er’s acetate decreased WB viscosity more than albumin, ACD-plasma, dex- 
tran 70 and gelatin in this order both on blood diluted in vitro and in vivo, 
Since the duration of volume expansion in vivo after Ringer’s solution was 
very transient the absolute decrease of viscosity was less than that obser- 
ved after all other solutions except gelatin. P viscosity increased with dex- 
tran 70, was unchanged with ACD-plasma, dextran 40, and gelatin and de- 
creased with albumin and Ringer’s acetate dilution, A marked red cell ag- 
gregation occurred with dextran 70 and gelatin with increasing concentra 
tion. No aggregation of red cells was seen with dextran 40, albumin, or 
Ringer’s acetate, It is concluded that lasting hemodilution in vivo is best 
achieved with albumin and dextran 40 based on their viscosity influences, 
volume expansion duration and their lack of red cell aggregation properties, 


INTRODUCTION 
Periferal vasoconstriction and cellular aggregation leading to venular stasis and microves- 
sel obstruction are pathophysiological effects seen in various kinds of shock, resulting in 
reduced and uneven blood flow through the microcirculation (1-7), Surgical trauma and 
shock is associated with loss of proteins and water to the extravascular space, resulting in 
decreased circulating plasma volume and elevated hematocrit (8, 9). Hemoconcentration is 
a common cause of an increased viscosity of blood, Surgical trauma induces red cell aggre- 
gation by protein changes, mainly by increased fibrinogen concentration with increased 
blood viscosity especially at low shear rates (10, 11). These factors cause a decrease of 
venous return and cardiac output, An uneven blood flow distribution causes hypoxia and tis- 
sue acidosis (7, 12-14), Hemodilution is a therapeutic method in shock and other conditions 
associated with hemoconcenfration, It will counteract the detrimental pathophysiologic de- 
fects in such conditions, 
The purpose of these experiments was to compare the effects of different plasma substitutes 
on blood viscosity and red cell aggregation when used for hemodilution and when given in 
comparable volumes and doses, 


MATERIAL AND METHODS 





In vivo experiments 





Sixty mongrel dogs of both sexes, weighing 20.0 + 3.8 kg, were anesthetized with keta 
initially 4 mg/kg, and continued with 1 mg/kg/hour, Pancuronum bromide was used a 
muscle relaxant at a dose of 1 mg/kg/hour, The animals were intubated and ventilated 

room air, Catheters were inserted through the right carotid artery into the aorta 

ing and through the external jugular vein for injection of fluids and drugs. After « 

surements shock was induced by laparotomy with splenectomy and exteriorization of 

small intestine and placed in wet gauze for three hours, Then new measurements were made, 
The intestine was thereafter returned to the abdominal cavity and the abdomen closed, The 
dogs were randomized into eight groups. The test fluids were warmed to 37 C and giv 
during a period of 20 minutes, The animals were followed for another four hour period, 


en 
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In vitro experiments 





Samples of human ACD-blood (hct = 43 %) were used for in vitro dilution experiments, Di- 
lution to hct 30 was achieved by addition of 3 ml of plasma substitute to 7 ml of the ACD 
whole blood, Erythrocyte sedimentation rate (ESR) was examined for samples of whole 
blood (hct = 30) with increasing concentrations of colloid. The measurements were made 
in plastic tubes with a diameter of 3 mm and read after 60 minutes, Red blood cell aggre- 
gation was assessed in a hemocytometer chamber. Whole blood was diluted with serum 
1:200) and the sample was added to the hemocytometer, Each monodispersed red cell or 
each red cell aggregate was counted after 10 minutes, at seven different concentrations of 
colloid from seven different blood donors, 


Measurements 





Whole blood (WB), plasma (P), and packed cell (PC) viscosities were estimated on 1 ml 
samples of heparinized blood diluted in vivo or on ACD-blood diluted in vitro at 37°C, 

using a Wells- Brookfield cone plate viscometer model LVT. After WB viscosity had been 
determined, plasma and red cells were separated by centrifugation at 2000 G for 15 minutes 
Viscosities were then determined on plasma and on packed cell mass separately. WB visco- 
sity was determined at shear rates 230, 115, 46,,23, and 11,5 sec™*, P viscosity at 115 
sec °, and PC viscosity at shear rate 11,5 sec °, Hematocrit was measured in duplicates 
in microhematocrit tubes after centrifugation at 12000 G for three minutes. No correction 
for trapped plasma was made, 


Test solution and experimental groups 

Group N: Eight animals were shocked, the intestine was returned to abdomen but no infu- 
sion was given, Group R: Eight animals were infused with Ringer’s acetate, (Content: Na* 
134 mEq/l, Ca“ 1.8 mEq/1,Cl 111 mEq/l, acetate 31 mEq/l, osmolality 260 mOsm/1, 
pH 7.0). Group G: Eight animals were infused with gelatin (polypeptides) with an average 
weight molecular weight (Mw) 5000 and an average number molecular weight (Mn) = 
6000, Group D: Eight animals were infused with dextran 40 of Mw = 40000 and Mn = 26000. 
Group M: Eight animals were infused with dextran 70 of Mw = 70000 and Mn = 41000. 
Group P: Eight animals were infused with dog ACD-plasma. Group A: Five animals were 
infused with albumin Cohn, made from dog plasma according to the Cohn fractionation me- 





thod, modified by Bjérling (15), (Content: Total protein 41 g/l, approximately 90 % of 
which was albumin (AB KABI, Stockholm, Sweden), Group E: Seven animals were infused 
with albumin PEG, made from dog plasma by fractionation with polyethylene glycol (PEG), 
according to the method of Bjérling (16). (Content: Total protein 41 g/l, approximately 

80 % of which was albumin (AB KABI, Stockholm, Sweden), 


All colloids were given as 3,5 % solutions in saline at a dose of 1.5 g/kg, corresponding 
to 43 ml/kg. Plasma was given in the same volume, Ringer’s acetate was given in a volume 
three times larger, or 129 ml/kg. 


Statistical methods 





Standard statistical methods were used to compute mean values and standard deviations, 

The effect of infusion is defined as the mean difference between pairs of individual values 
at shock and values at different times after reposition of the intestine. The comparisons 

of different groups were carried out as univariate analysis of variance with the shock ef- 
fect as a co-variate (17, 18). 


RESULTS 


The viscosity of different plasma substitutes dissolved in saline increased with increasing 
concentration (Fig. 1). Normal heparinized plasma had a viscosity of approximately 1,3 - 
1,4 cp, ACD-plasma had about 15 % lower viscosity or 1,2 cp because of its dilution with 
the ACD-solution, At a colloid concentration of 3.5%, albumin had a viscosity of 0,8 cp, 
dextran 40, and gelatin 1.3 cp, and dextran 70 1,6 cp, Gelatin in concentrations higher than 
% was not available, In these experiments saline had a viscosity of 0,7 cp. 


; » 


Plasma viscosity (P) 
otal plasma protein concentration and albumin concentration decreased slightly to 96 % 


and 91 %, respectively, of the initial values (Table 1), Plasma protein and albumin concen- 
trations changed after infusion in relation to the hemodilution effects of the various fluids, 


P viscosity did not change during shock (Table 2), P viscosity increased after infusion of 
dextran 70 (p <0,001), Ringer’s acetate lowered P viscosity (p< 0,001). These influences 
were most pronounced one hour after infusion but were still effective at four hours after in- 
fusion (p< 0.05), Gelatin, dextran 40, and ACD-plasma did not change P viscosity signifi- 
cantly, while albumin Cohn caused a slight decrease, Similar relative changes were ob- 
served on samples diluted in vitro (Fig. 2). 
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FIG, 1 FIG. 2 


Apparent viscosity versus concentration Change in plasma viscosity after dilution 
for different colloid solutions. in vitro and in vivo with 3,5 % colloid and 
Ringer’s acetate solutions, Mean + SEM, 


Legends as in Fig. 4. 


TABLE 1 


Total Plasma Protein, Albumin, and Hematocrit Changes during Shock 








—— : : 
Control Shock Effect of shock 


NS ea 


——EEE 








Total protein g/l -64+6,9° 60,3 + 6. 2.4 + 5.6 * 
Albumin (g/1) + 2.3 23.0 ‘ -2.3 .i ** 


Hematocrit (%) +5.,7 790 8 + De 6.2 + 4,2 * 























* Mean + SD. 


Whole blood viscosity (WB) 


WB viscosity increased during shock at all shear rates, more so at lower shear rates 
(Fig. 3 and Table 2), concomitant with the increase of hematocrit from 50,9 to 59.1 %, 
All infusions decreased WB viscosity because of hemodilution. When comparing WB vis- 
cosity change at equal hematocrit decrements, Ringer’s acetate lowered WB viscosity 
more than albumin Cohn, ACD-plasma, dextran 40, albumin PEG, dextran 70, and gela- 
tin in this order, When blood was diluted in vitro with 3,5 % solutions the same relative 


influences were obtained with the different plasma substitutes (Fig. 4). 





Packed cell viscosity (PC) 
PC viscosity increased slightly during shock to 106 % of control values (Table 2), and was 
further elevated with dextran 70 and gelatin infusions in relation to all other groups 


FIG. 3 


entipose 


Apparent whole blood viscosity before and 


after three hours of intestinal shock, 
Mean + SEM. 
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p <0,.001). PC viscosity was unchanged after dextran 40 and ACD-plasma and decreased 
with albumin or Ringer’s infusions. PC viscosity changes were similar in blood when di- 


luted in vitro with colloids of 3,5 concentrations (Fig. 5, Table 2). 
TABLE 2 


Viscosity Changes during three hours of Intestinal Shock 








Shear rate, sec | Control S Effect of shock 
aden cinenameeieniat EEE 





230 é ,° ; + 1,05 + 0,87 ** 


115 i - + 1,02 + 0,93 ** 


WB Viscosity ‘ - . u 


+ 1,20 
+ 1,67 4 





P Viscosity 





PC Viscosity 

















* Mean + SD 


Change in whole blood viscosity in relation 

to hematocrit change after dilution in vitro 

and in vivo with 3,5 % colloid and Ringer’s 
ate solutions, Mean + SEM, 


Erythrocyte sedimentation rate (ESR) in vitro 
ESR increased with increasing concentrations of dextran 70 or gelatin, and was unchanged 
with increasing concentrations of dextran 40 or albumin (Fig. 6). 
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FIG. 5 FIG. 
Change in packed cell viscosity after dilu Erythrocyte sedimentation rate 
tion in vitro and in vivo with 3.5 % colloid versus concentration for different 
and Ringer’s acetate solutions, Mean + SEM, colloid solutions. 


Legends as in Fig. 4. 


Microscopic red blood cell aggregation count in vitro 


Approximately 50 % of red blood cells were aggregated in rouleaux formations in heparin 
ized plasma and 20 % in serum, With increasing concentrations of dextran 70 or gelatin the 
number of free monodispersed cells decreased and the number of aggregates increased up 


to 12 gZ colloid substance per liter above which concentration this system seemed to reach 
’ y 
Figs. 7 and Oo}. 


a saturation point, where no further cells were available for aggregation 


FIG, 7 
Monodispersed red blood cell 
DEXTRAN 4 count versus increasing colloid 
concentration. Mean + SEM, 
S serum; P = plasma, 


q 12 5 


CONCENTRATION 
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Albumin caused no aggregation, Dextran 40 did not aggregate at concentrations below 
i2 g/l. Exceeding this concentration a slight aggregation was seen but less than with 
normal plasma, 


——* GELATIN 
} DEXTRAN 70 


7’ DEXTRAN 40 


+ 
~>+——-?_ ALBUMIN 


g 12 15 18 
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FIG, 8 
Red blood cell aggregate count versus increasing 


colloid concentration. Mean + SEM. S serum; 
P = plasma. 


DISCUSSION 


The purpose of the present experiments was to study blood viscosity and red blood cell 
aggregation at comparable hemodilution and concentrations of colloids. Colloids have dif- 
ferent inherent viscosity and aggregation effects depending upon their physicochemical 
properties (Figs. 1, 6-8). The in vitro viscosity of plasma and of whole blood after dilu- 
tion will therefore be determined in a predicted way by the concentration and volume of a 
specific substitute used for dilution. Such in vitro experiments do not apply to the in vivo 
situation, where the hydration status further determines whether water will move in or 
out of vascular space, depending on the hydration of the recipient and the colloid osmotic 
activity or concentration of the fluid given (9, 13). 


The different plasma substitutes had similar relative influences on blood viscosity both 
when diluted in vivo and in vitro, when compared at equal hematocrit change and concen- 
tration of colloids. The greater general decrease in whole blood viscosity in vivo per 
change of hematocrit was explained by a viscosity decrease even when no infusion was 
given (group N) (Fig. 4). 


During shock total protein concentration did not change substantially which is consistent 
with unchanged plasma viscosity during shock, 


Viscosity measurements in vitro reflect only a few flow determinants that are effective in 
vivo, Even if one accepts in vitro viscometry data as indicative of a trend, the pricise in- 
fluence is difficult to analyse in quantitative terms, The peripheral vasoconstriction and 
cellular aggregation with blockage of small vessels resulting in uneven or maldistributed 
blood flow is probably of far greater importance and is only to a minor extent determined 
by in vitro viscosity measurements (1, 6, 7, 12). 


Whole blood viscosity at higher shear rates (above 115 sec”) is mainly affected by hema- 
tocrit and the red cell deformability. At lower shear rates (almost zero) red blood cell 

aggregation plays an increasing role and causes blood viscosity to increase (19, 20), Du- 
ring shock, flow rates near or at zero exist, especially on the venous side of the capilla- 
ry. This stasis causes increased viscosity and increased resistance to flow, resulting in 


elevated capillary pressure and efflux of fluid and proteins, with further loss of plasma 
volume (21), 


This effect is further potentiated by hemoconcentration seen in shock, 


Packed cell viscosity is usually used in order to eliminate hematocrit variations and is be- 
lieved to reflect direct cell to cell interactions and red cell deformability (20), The packed 
cell viscosity changes in the present study were similar when blood was diluted in vivo 
and in vitro. In addition, the changes in packed cell viscosity were consistent with those 
seen in erythrocyte sedimentation rate and in the microscopic aggregation test for the dif- 
ferent plasma substitutes, suggesting that packed cell viscosity is influenced by changes 

in red cell aggregation, 


Red blood cell aggregation occurs regularly after injury (4-6) and contributes to elevation 
of whole blood viscosity at low shear rates (19, 20). The degree of aggregation is directly 
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related to protein changes, mainly fibrinogen and globulins (22). Increased aggregation is 
reflected and measured by increased erythrocyte sedimentation rate (23), The in vivo con- 
sequence of cellular aggregation is capillary blockage, maldistribution of blood flow, tis- 
sue anoxia and tissue acidosis (1, 6, 12-14), Macromolecules, other than proteins, also 
cause red blood cells to aggregate. The difference between dextran 70 and dextran 40 is 
striking. The aggregating capacity depends on the concentration, the average molecular 
weight and the molecular weight distribution (23-26). The results in the present study re 
garding viscosity and red blood cell aggregation are consistent with those reported by 
Brooks, Jan, Chien, and Lim (24-27). 


CONCLUSIONS 





At dilution to equal hematocrit Ringer’s acetate decreased whole blood viscosity more 
than albumin Cohn, ACD-plasma, dextran 40, albumin PEG, dextran 70, and gelatin, 
in this order. 


Albumin Cohn and Ringer’s acetate decreased plasma viscosity. ACD-plasma, dextran 
40, and gelatin caused no change, and dextran 70 increased plasma viscosity. 


These relative influences on blood viscosities from plasma substitutes were similar 
when blood was diluted in vitro and in vivo. 


Albumin, dextran 40 and Ringer’s acetate did not cause red cell aggregation in con- 
centrations less than 12 g/l serum, Dextran 70 and gelatin caused marked red cell 
aggregation at corresponding concentrations. 


It is concluded that lasting hemodilution in vivo is best achieved with albumin and dex- 


tran 40 based on their viscosity influences, volume expansion duration, and their lack 
of red cell aggregation properties. 
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A COUETTE VISCOMETER FOR SHORT TIME SHEARING OF BLOOD 
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ABSTRACT 


A Couette viscometer with axial flow is proposed which allows blood to be exposed to 
defined shear stresses of short duration. A prototype of this viscometer is tested. By 
measurements and simple calculations the range of application is established. Measurements 
with blood show the possibility to control the flow conditions which may lead to haemolysis. 


INTRODUCTION 


The damage of blood particles through shear stresses can be of essential importance for artificial organs. 
Not only the magnitude of the shearstresses but also the time they act on the particles is important. 
Various experimental set up were constructed through which both the influence of magnitude and time 
can be investigated. Figure 1 shows a summary of the ranges examined. Williams (1) and Champion et al. 
2) use a cone and plate system. Leverett et al. (3), Nevaril et al. (4) and Whiffen et al. (5) carried out 
measurements in a combination of a cone-plate and a Couette system, Bacher et al. (6) and Keshavia et 
al. (7) in capillaries. Lambert (8) examined the deformation of erythrocytes in a narrow slit. Short loading 
times were reached by Rooney (9) with oscillating gas bubbles, by Williams et al. (10) with oscillating 
wires, by Bernstein et al. (11), Blackshear (12) and Forstrom et al. (13) with turbulent jets. Short 
shearing times are of particular interest, since in natural and artificial circulations hardly a constant 
shearstress will develop for minutes as in (1,2,3,4,5). Naumann (14) arrived at very short loading times of 
about 10 ~ sec for the flow through an artificial heart valve or through a stenosis. 
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Fig. 1: Range of shear stresses and shearing times examined and 
test devices used by other authors 
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Well defined flow conditions are required for the investigation of shear induced haemolysis. The shear 
stress distribution as well as the duration of the load should be known. This is not the case in many 
experiments (9,10,11,12,13) either the estimate of the shear stress is not sufficiently accurate or the 
particles may be loaded repeatedly in the vortex motion of the flow. Clearly defined conditions are found 
in capillaries (6,7), but here the shear stress is not constant. In a slit (8) only one erythrocyte can be 
observed. The existing variation in its volume is taken as a measure of the amount of penetrated 
haemoglobin. 


In this paper a Couette system is described, in which high shear stresses can be generated in short loading 
times. Both parameters can be varied independently. For this purpose the Couette system has an axial 
flow. The rotation of the inner cylinder controls the shear stress and the axial velocity regulates the 
loading time. This system is described and tested, and first measurements with blood will be presented in 
the following. 


DEVELOPMENT OF A COUETTE SYSTEM WITH AXIAL FLOW 





Reproducable shear stresses can be generated in the Couette system by the choice of the speed of 
rotation of the rotating inner cylinder, if the geometry and the viscosity of the medium are known. The 
diameter of the inner cylinder d: 58mm and a gap width b=0,115 mm produces a shear velocity 
of # = 2.1+10° sec © for its maximum rotational speed n = 8000 R/min. The height of the cylinder is 
h = 5mm. A proportioning pump produces a well known volume flow V. Since the geometry is known, the 
shearing time t,, results to t h-it-d.-b/V. A sectional drawing of the system is shown in figure 2. Blood 
flows through tRe inlet A tangential into a toroidal channel with stationary walls, from here through the 
shearing region of the Couette system and from there into another toroidal channel with stationary walls 
and finally through the outlet B, where the probe is taken. The whole device is thermo-requlated at the 
outer cylinder and within the apparatus. The inner cylinder is fastened at the shaft which can be adjusted 
in axial and radial direction. 


cooling 
= 





inflow A 
—— 


Fig. 2: Sectional drawing of the Couette visco- 
meter 


Since high shear stresses should be developed only within the Couette system by the rotating inner 
cylinder, this cylinder is stuffed against the stationary hgysing (sliding velocity V@ 424 m/sec). A 
suitable packing was found in an axial wheel shaft packing. “ The inner cylinder, consisting of backing 
coal, reinforced by an outside shrinked on stainless steel (VA) ring, moves on standing coale rings. The 
behaviour of the packing, the abrasion of which is indicat@d with 1 mm per 10 working hours will be 
described. The toroidal channels and the in- and outflow tubes were made out of titanium. 


The flow in a Couette system with axial flow has been investigated by many authors. Chung and Astill 
15) theoretically and Gravas and Martin (16) examined experimentally the onset of Taylor vortices 
depending on the mean axial velocity ¥ and the rotations per unit time n of the inner cylinder. The flow 


*) This packing was delivered by Fa. GOETZE WERKE, 567 Opladen 
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modes which developed are shown in figure 3. 


turbulent flow 


_— 
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velocity 


turbulent flow 


fl 
— a with Taylor - vortices 


laminor flow 
with Taylor -vortices 





Oo mean oxial 


rotations per unit time n 
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Fig. 3: Modes of flow which may develop 
in a Couette system with axial flow 


In the system proposed here, neither the axial nor the tangential velocity profiles are fully developed. 
Astill et al. (17) investigated this flow experimentally and Coney and El-Shaarawi (18) theoretically. 
Velocity- and pressure distributions were obtained for the developing flow region, however not in a range 
of axial velocities and rotational speeds which are of interest here. Estimate of the two limiting cases, 
either rotating flow without axial flow or axial flow only show, that both profiles are developed at least 
after one tenth of the length which is given by the construction. Two requirements should be fulfilled: a) 
shear stresses caused by the rotation have to exceed those caused by the axial flow, b) haemolysis 
measurements should be done under laminar flow conditions only. These conditions establish the range, in 
which the system can work meaningful. This range is shaded in figure 3. The shear stresses given in some 
of the following figures are those caused by the rotation only and are thereby to be considered as the 
lower limiting case. 


Since the axial flow causes a velocity profile, the residence time of blood particles in the gap will depend 
on their distance from the cylinder. An erythrocyte which is in a distance of 104sm from one cylinder 
will remain twice as long in the slit as one which is in the middle. The loading time ty is the mean 
passage time through the gap. 


First haemolysis measurements were done in the loading range which is shaded in figure 1. This range can 


be extended by changing the geometry of the gap by the interchangeable outer cylinder or by varying the 
viscosity of the test fluid. 


TESTING THE SYSTEM 





The Couette viscometer was tested repeatedly in order to ensure a reproducable and well defined flow 
field. The influence of possible error sources were estimated. For the operation of a Couette system it is 
essential that the cylinders are concentric. As was shown by Ullrich (19), the pressure loss caused by 
axial flow depends on the excentricity e of the cylinders (see figure 4). Measurements of the pressure loss 
show a good agreement between theory and experiment. This method allows only a rough checking as to 
whether or not the cylinders are concentric, since the differences in the measured pressure losses at 
e = 0 are small. A very sensitive method is to measure the static pressure at three locations uniformly 
distributed on the perimeter. Excentrically rotating cylinders will give an asymmetrical pressure 
distribution. Without axial flow the rotating inner cylinder is shifted relative to the outer cylinder and 
will provide a symmetric pressure distribution. 


In- and outflow channels are symmetric to a plane normal to the axis of rotation. Different pressures in 
these channels at rotating inner cylinder without axial flow indicate an asymmetry to this plane. The 
inner cylinder is too high to too low with respect to the outer cylinder. The results of the measured 
pressure difference 4p between in- and outflow channel depending on the axial shifting distance s of the 
inner cylinder are shown in figure 5. The curves for the two directions of rotation are caused by the 
backlash of the bearing. After fixing the sense of rotation (inflow in direction of rotation) the pressure 
difference Ap _ is adjusted to zero. 


Laminar flow should develop in the whole range of operation. In this case the axial pressure difference 
between in- and putflow channel depends linearly on the loading time t,,. As figure 6 shows, this is true 
for #@15.53 +10 sec ~. It is concluded that for the pressure loss in the gap quadratic terms may be 


neglected. The deviation from this behaviour at higher shear velocities can be explained by a small 
remaining pressure difference even if the axial flow is zero (see figure 5). The straight lines are shifted 
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Fig. 4: Pressure loss of axial flow as a function 
of the excentricity e of the cylinders 
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Fig. 5: Test of the axial symmetry. Clockwise or anticlockwise 
rotation when viewed from above 


with increasing shear velocity at constant loading time to higher pressure differences. 


As was shown by Karle (20), temperature influences the rate of haemolysis, too high temperatures 
jamage the erythrocytes directly. Therefore the temperatures occuring in the Couette system should be 
known. The temperature is measured by means of a probe introduced through a pressure tap next to the 
fluid penetrating the gap. The rise in temperature is plotted against the time and the stationary end 
temperature is transferred into a r- t..-diagram. Points with constant temperature increase 4 Pare 
onnected and shown in figure 7. These Curves are obtained without additional cooling in the thermostat 
reservoir. The whole system is heated up by the frictional heat produced by the axial wheel shaft 
packing. With additional cooling the ternperature rise can be halfed. It was found to be independent of 
the viscosity 7 of the test fluids (1 cP << 10 cP). The measurements with blood were carried out with a 
temperature of 2?= 30° C in the reservoir and additional cooling. 


Because of the high shear velocities which are required, the angular velocity of the inner cylinder has to 
be increased to 800 sec ©. In this case a pressure gradient within the cooling fluid inside the housing of 
the system is to be expected which may cause a leakage of cooling fluid through the axial wheel shaft 
packing into the blood. The leaking water volume was measured as a function of the shear velocity. 
Especially for measurements with blood not the absolute value of the leaking volume seems to be of 
interest but rather the percentage of volume leaving the outflow B coming from the cool jng circuit. This 
; plotted in figure,8 as leakage in percent against the loading time. F or *<4-10° sec © the packing in 
tight, for 725-10 sec ~ the leakage volume is increased. The maximum loading time was 0.69 sec for 
experiments with blood which results in a leakage of 7 percent. To avoid osmotical haemolysis 
physiological solution of sodium chloride is used as cooling fluid in experiments with blood. Experiments 
with different coloured cooling water show the leakage to be apportioned half by half on each of the 
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Fig. 6: Pressure loss in the gap by axial flow and rotating Fig. 7: Increase in temperature in the gap 
inner cylinder without additionalcooling 
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Fig. : Leakage throug . al wheel shaft packing 


In experiments with blood it is necessary to rinse the blood leading system with physiological solution of 
sodium chloride and to keep the system filled with it to avoid blood air contact. The sodium chloride 
solution in the system is to be displaced by blood conveyed by the proportioning pump. The time required 
for this process is to be determined as a function of the operating conditions. Instead of the sodium 
chloride solution the system is filled with coloured water. The rinsing process is determined by measuring 
the colour goncentration with a photometer. The process is completed if the concentration of the colour 
is below 1 ~/o. The time which is known now allows one to calculate the volume to be displaced to get 
blood with known composition and load at the outflow. This volume V is plotted as a function of the 
loading time t,, in figure 9. Differences between axial flow with and without rotation exist only for 
t..> 0.1 sec. Thé displaced volumes are independent from the shear velocity for axial flow with rotation. 
us displaced volumes cannot be used in experiments with blood. 


The possible range of operation is demonstrated in figure 10. This range is limited, by the maximum 
rotational speed which corresponds to a maximum shear velocity of = 21.37 +10 sec. Since the 
maximum wall shear stress 7 should not increase above 10% of the shear stress caused by the 
rotation T the possible loading time t,, is limited. Points of measurement on,the left side of the 
curve T oi? are to be excluded. As an example a point x 15-10 sec andt 0.1 sec is 
considefed. The diagram shows: a) the leakage of cooling water into the flowing blood is smaller than 
1%, b) the temperature incrase AYwill be 7 C without additional cooling, with cooling an increase of 


rat’) 


3.5° C is to be expected, c) Taylor vortices will not develop in fluids which have a viscosity of > 2.3 cP. 
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Fig. 9: Displaced volurne to get blood with known composition 
and load 
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Fig. 10: Possible range of operation 
ynstant increase in termperature (see Fig. 7) 
ynstant leakage (see Fig. 8) 
mset of Taylor vortices depending on the vis- 
sity 7 of the fluid (Astill (19)) 


HAE MOL YSIS MEASUREMENTS 





The measurements of blood damage described below were done with fresh heparinised pig blood. The rate 
of haemolysis Hg is defined as the ratio of plasmahaemoglobin released during the load 4 Hb to the total 
amount of haemoglobin Hb in percent. The plasmahaemoglobin concentration at the beginning is PHb,and 
at the end PHb of the experiment, the haematokrit Ht and the whole blood haemoglobin concentration 
VHb are to be measured at every sample. These values are evaluated with the formula 





Ai - PHb - PHb 
A AM ry. (1-Ht/100)( o) 


Hg = - 700 


Hb VH6b 


The blood damage thus defined is plotted in figure 11 depending on shear velocity and loading time. In a 
pilot experiment could be shown that there was no blood damage detectable with axial flow alone. All 
curves show the same shape at constant _—- fime, although t,, is varied over more than two orders of 
magnitude. For t,<3.4-10 “sec and j<3- sec - no haemolysis is detectable. The shape of the curves is 
similar to those obtained by Williams (1) who varied the shear stress at constant loading time of 5 
minutes. 


A diagram of haemolysis, in which the shear stress is plotted depending on the shearing time with the 
rate of haemolysis as parameter, can be obtained from figure 11. The viscosity of blood was chosen to 
be 7 = 3 cP. The diagram of haemolysis (figure 12) illustrates the influence of time on haemolysis. F or 
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Rate of haemolysis depending on shear Fig. 12: Diagram of haemolysis 
velocity and loading time. Haemolysis 

for the operation of the loading pump 

alone is not detectable. 


one shear stress a great intervall of time is needed to reach, beginning with Hg = 0, total haemolysis. 
Rates of haespplysis Hg>5% are represented by straight lines which are proportional to the loading 
quantity T-t as was already shown by Lambert (8). This dependence also results from a theoretical 


investigation by Richardson (21). 


SUMMAR Y 


A Couette viscometer with axial flow to generate high and clearly defined shear stresses and short 
shearing times is tested, for the projected range of operation. The concentricity of the cylinders is 
controlled by measuring the static pressure at three locations uniformly distributed on the perimeter. 
The symmetry to a plane normal to the axis of rotation is tested by the measurement of the static 
pressures in the toroidal channels with the inner cylinder rotating. With additional axial flow the results 
of the measurement of the pressure difference in these channels permit to neglect quadratic terms in the 
pressure loss in the gap. For the measurements with blood the increase in temperature, the leakage 
through the packing and the displaced volume which is required to get blood of known composition and 
load at the outflow are determined. The results of these measurements are summarized in the operating 
diagram. Experiments with pig blood show the applicability f the viscometer to measurements of blood 
damage by shear stresses. Shear stresses T<7-~-10° dyn/cm” and shearing times th? 3-10 ~ sec can be 
obtained. 
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ABSTRACT A capillary viscosmeter is described which allows measurement of 
the shear rate dependency of fresh blood viscosity at low shear 
rates. The instrument has been used in a preliminary study to 
evaluate differences between fresh and EDTA anticoagulated blood. 
In a shear rate range of 0.5 to 1.3 sec l the viscosity of fresh 
blood is more strongly dependent on hematocrit than is the viscos 
ity of EDTA anticoagulated blood. Also, the viscosity of anticoag 
ulated blood is more sensitive to day-to-day variations than is 
that of fresh blood. For some subjects, such variations may amount 
to as much as 30%. 


INTRODUCTION 


There is a considerable body of literature on the rheological proper 
ties of blood, both normal and abnormal. At low shear rates whole blood i 
pseudoplastic fluid. Furthermore, it has been well documented that both the 
absolute value and shear rate dependency of whole blood viscosity vary sig 
nificantly in various abnormal or disease states (1). 

Almost the entire literature dealing with whole blood rheology is based 
upon studies performed on anticoagulated samples. A number of investigations 
using fresh blood have given conflicting results: several conclude that there 
is no effect of anticoagulants on blood viscosity, while others suggest that 
there are differences in the viscosity of fresh vs. anticoagulated blood 
(2-5). Because most precise viscometric techniques require considerable set 
1 to one shear 


7 
Thi 


S a 


up time, fresh blood viscosity measurements are usually limite: 
rate value for each sample before significant clotting takes place 
paper describes a capillary ramp-response viscometer which c: determine 
whole blood viscosity as a function of shear rate on a ingle ample of fre 
lood. feasurements over a wide shear rate range ci ye complete: 1 le 
blood Measurement ver a l I t ran an | ympleted in 1] 
than one minute after venipuncture. The method is capable of giving shear 
stress/shear rate data directly, giving true viscosity, as well as apparent 
viscosity, and is not dependent upon the a priori assumption of a specific 
rheological model. 


DESCRIPTION OF THE INSTRUMENT 


A block diagram of the instrument is given in Figure l. The unusual 
feature of this capillary viscometer is that it is operated 
flow device, i.e., the flow rate is deliberately varied 
tion of time. the version described here the flow 
function and is in the pseudo-steady state regime; that i 
large compared to the time constant of the system. The flow ramp is 
by using a modified syringe pump (Harvard Model 990 Servo-Control Sy 
Pump). A linear voltage ramp is generated by a suitable function gene 
This signal drives the syringe pump controller, which then generates a 
output linearly varying with time. The ramp function generator output 
drives the arm of an X-Y recorder, so that the X axis records a signal 
rectly proportional to flow rate. The pressure drop across the capillary is 
measured by a differential transducer (Validyne Model GP-45 Transducer with 
Model CD-12 Transducer Indicator). The output of the transducer is re i 
as the pen signal on the X-Y recorder. When a Newtonian fluid flows throug! 
the capillary, the result should be a linear plot whose slope is proportional 
to the viscosity of the fluid. When a non-Newtonian fluid flows through the 
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Block Diagram of Capillary Ramp Function Viscometer 


capillary, a nonlinear plot should be obtained whose slope at any point is 

proportional to the apparent viscosity at the average shear rate character- 
ized by that particular flow. Ising the Rabinowitch equation (6), the data 
at different flow rates can be transformed to give an absolute viscosity. 

Characteristics of the instrument were studied and calibration was 
carried out using several sucrose solutions of varying viscosities (5 - 40 CP) 
as well as anticoagulated human blood. The pressure transducer was calibrated 
against a micromanometer. The viscosity determined from the ramp function 
response agreed well with values determined from steady state runs, indicating 
that the pseudo-steady state assumption was valid over the range of flows and 
viscosities used. Analysis of the transient characteristics of the system 
indicated that with the capillary and transducer used the time constant of 
the system was about 1 second for a viscosity of 24 CP. Experiments indica- 
ted that the optimum ramp times were onthe order of 10 seconds; this gave a 
linear response over most of the output curve (for Newtonian samples) and 
failed to demonstrate any changes with anticoagulated blood that could be 
attributed to sedimentation or aggregation. Further details of the instru- 
mentation and calibration can be found in Downs (7). This instrument is a 
modification for low shear rat of that developed by Kron (8). 

To carry out a test, a constant initial flow was set and a decreasing 
ramp of flow vs. time was run. Data from a calibration run of this type are 
shown in Figure 2, for a sucrose solution. Overall time of the experiment 
is 10 seconds, starting from an initial shear rate of 5 sec™+. There is an 
initial curvature of the plot due to the compliance of the pressure trans- 
ducer. The curve then becomes linear, as expected of a Newtonian fluid. 

When the impressed flow falls to zero, however, the pressure is still above 
zero. This is due to the displacement of the transducer membrane; because of 
the high viscosity of the fluid in the capillary, it takes a finite time for 
the pressure to decay to zero. This behavior can be understood in terms of 
the following analysis. For Newtonian fluid, assuming the Hagen-Poiseuille 
law for flow in a capillary obtains, we have 


8uLQ - ; 
, ° 5 (1) 


P = 
Ro 
where P = pressure difference across the capillary of length L and radius R: 
and yw is the dynamic viscosity. The flow rate through the capillary Q has 
two components, as illustrated in Figure 3. Qp is the induced flow due to 
the syringe pump; Qt is the flow due to the displacement of the transducer 
diaphragm. The components are given by 
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i 
‘pl 
-dV 
Q = t 
t dt 


where Vt is the volume of the transducer. The total flow through the capil- 
lary is Q = Qp + Qt. It is assumed that the sidearms A and B are designed 
so that pressure drops in themare negligible. The flow due to the transducer 
is linearly dependent on the pressure measured by the transducer through the 
compliance k]: 
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the Ramp Response of a Newtonian Fluid 
(24CP sucrose solution) 


Combining these equations, we have 


k,{Qpj(1 - t/te) -k, o) 


te being the time at which the pump is turned off. Solving this equation 
gives, for the pressure vs. time: 


From Equation (1), Po, the initial pressure before starting the ramp, is 
given by k2uQpi- Using this, we can solve Equation 6 explicitly for the 
viscosity, giving 


, P/P, . t/t, 2 } 
‘J - exp (-t7k) ku) 
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For the parameters used the present apparatus, after a few seconds the 
r r , 
exponential term decays yivil 


e ‘nit ; > 
kjk, f _ t/t. (8) 


is the Newtonian viscosity (or apparent viscosity for non-Newtonian 
This explains the linear portion of the curve in Figure 2. Further- 
is apparent from Equati (8) that at t = t,, i.e., at the end of the 
pressure is not equal zero, but is given by a finite value which 


(9) 


From Equations (7) and (9), it is seen that the term Ky kau is the first order 
time constant of the system, and that the nonlinear effects can be minimized 
by making the time constant small compared to the ramp time. Within the con- 
straints of maintaining transducer response and low shear rate, it was pos- 
sible to obtain useful linear data over about 80% of the ramp. It should be 
noted that since kj is the onl compliance term used, it effectively includes 
the entire system compliance. 

For a non-Newtonian fluid, one can apply a similar analysis except that 
the middle portion of the plot will now also be curved because yp = £(Q). At 
each instant of time, an apparent Newtonian viscosity Uapp can be calculated 
from Equation (8) which assumes the applicability of Poiseuille's law; the 
corresponding flow rate is obtained from Equation (2). From this and the 
pressure one can then calculate the Newtonian wall shear rate Ywn and the 


wall shear stress Tw: 


YWN 


Using these data and the Rabinowitch equation, one can obtain the true 
non-Newtonian viscosity any at any shear rate via: 
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Figure 4 shows a raw data plot fer whole blood. The slight oscillations 
in the curves are due to residual noise from the syringe pump. 

For the studies carried out here, a capillary of 2.108 mm diameter and 
24.3 cm long was used. The usable shear rate range was from about 1.5 down 
to 0.1 sec-4. The shear rate range can easily be extended by varying the 
gear setting of the pump and is ultimately limited by the sensitivity of the 
pressure transducer. The temperature was maintained at 37° by immersing the 
entire capillary-transducer assembly in a water bath. 





PRESSURE 











Tracing of Actual Pressure-Flow Curves from X-Y Plotter. 
Fresh Blood Sample from Normal Subject. Curves 1-3 run 
on same sample approximately 20 seconds apart. 


EXPERIMENTAL STUDY 
Populations 


An experiment to test the utility of the instrument was designed. The 
apparent viscosity of fresh vs. anticoagulated blood at low shear rates was 
determined. Subjects were students ranging in age from 18 to 23 years; 
there were 15 males and 4 females. Ramp tests were run with shear rates 
ranging from 1.5 down t¢ 0.1 sec”~+- The data are reported at two shear 
rates, 1.3 and 0.5 sec™+. 

Because the results of previous work had indicated that significant 
viscosity variations could occur due to time of sampling and behavior of the 
subjects (time since last meal, stress from examinations, prior weekend par 
tying, etc.), the study was carried out under as controlled conditions for 
the subjects as possible. Data from each subject were collected on each of 
the three successive days. These were Tuesday, Wednesday and Thursday morn 
ings, which were selected as being the least stressful days of the students' 
week. A questionnaire was used to determine any relevant personal character- 
istics of the subject. Each subject was asked to come to the laboratory 
prior to breakfast and to refrain from eating or drinking anything after 
midnight of the night before. 
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Testing 


As indicated in Fig. samples can be introduced directly from the 
vein into the instrument, from a syringe. In the present studies, all 
samples (10 ml) were drawn from the antecubital vein to a syringe and then 
loaded into the instrument, since repeated trials indicated that the instru- 
ment could be loaded just as rapidly in this fashion. Prior to use, the capil 
lary was siliconized and the instrument flushed with blood to equilibrate 
all surfaces. It was then filled with saline. Experiments with anticoagu- 
lated blood demonstrated that about 4 ml of blood were needed to replace 
the saline in the capillary (volume 0.85 ml) to avoid dilution errors. 
Depending upon the skill of the operator, the instrument was ready for the 
initial run 20-30 seconds following withdrawal of the sample. Each run 
took 15 seconds, with 5 seconds between tests. Thus, in most cases, three 
or more complete replicate te could be carried out on the same sample 
in a period of less than 2 minutes after venipuncture. After loading the 
machine with fresh blood, the remaining sample was introduced into a vacu- 
tainer (of appropriate size for the blood volume) which contained the proper 
amount of dry EDTA. The anti gulated samples were run as soon as possible 
after the fresh ones, usually within one hour. Hematocrit was determined 
on the EDTA sample using the microhematocrit capillary tube method, on 
the same day. Between samples, the system was flushed with saline until 
clear. Also, the system was periodically cleaned more thoroughly by removing 
the capillary and reaming wit! large pipe cleaner, followed by recoating. 
Recalibration after extended use indicated negligible changes. 

In all, 57 blood samples drawn from 19 subjects were tested (one per 
day for each over 3 days). A mited number of tests were also carried out 
with samples from subjects showing various forms of vascular pathology, 
including stroke, TIA (transient ischemic attack), and sickle cell anemia. 
These tests were all made on EDTA anticoagulated blood, since it was not 
possible to transport the patients to the laboratory. 
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RESULTS 


A typical plot from a subject having a hematocrit of 50, in which three 
tests were carried out in succession, is shown in Figure 4. The tests are 
labelled 1-3 in the temporal order in which they were carried out. In this 
case the first test was completed in slightly under 1 minute. Each subsequent 
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test was completed 20 seconds later. There is a definite increase in the 
initial pressure Pp for each repeat, which are all carried out at the same 
initial flow rate O03. Therefore, there must be significant increase in 
viscosity over this time period, indicating that rheological changes associ- 
ated with coagulation may be taking place from the moment of venipuncture. At 
1 sec-4,the results indicate an 8% increase in viscosity in 20 seconds. Over 
the same time period, similar repetitions using the EDTA anticoagulated blood 
were reproducible, indicating that sedimentation was not causing the varia- 
tion seen in Figure 4. Allowing anticoagulation blood to remain undisturbed 
in the capillary for a long time (> 5 min) did produce changes that could be 
attributed to sedimentation. 
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Viscosity of Normal EDTA Anticoagujated Blood as 
a Function of Hemocrit at 0.5 sec” Solid Line: 
Least Squares Regression. Dashed Lines: 

95% Confidence Interval 


Because of this, in the present study, the earliest results were used 
as representative of the "fresh" blood sample viscosity. Actually,a better 
method might have been.to extrapolate viscosity values back to the moment 
of venipuncture, but this would have required much better control over and 
recording of samplestimes, with presumably individualized corrections for 
variations in coagulation rates between samples. Certainly further studies 
to quantitate this effect are in order. 

A semi - logarithmi¢ plot of fresh blood viscosity vs. hematocrit at a 
shear rate of 0.5 sec-+ is shown in Figure 5. A least squares regression 
line has been fit to the data. Also included are the 95% confidence limit 
bands for the individual measured values. A similar plot for the same 
samples after EDTA addition is given in Figure 6. The_results of statistical 
analysis of*these data and at a shear rate of 1.3 sec”™+* are given in Table l. 
The correlation coefficients for the regressions range from about 70 to 85%, 
indicating some scatter but highly significant correlations. These results 
indicate that the predicted slopes of the EDTA and fresh blood viscosity 
plots were significantly different at the 95% level at a shear rate of 1.3 
sec>+. At the lower shear rate, the difference in slopes was significant 
only at the 84% level. Interestingly enough, the two plots cross at hemato- 
crit around 45, indicating that measurements on hematocrits in that range 
would tend to show no differences between fresh and anticoagulated blood. 
However, at conditions of extreme high and low hematocrit, which might be 
typical of abnormal states, there appears to be a significant difference in 
the viscosity/hematocrit behavior of fresh vs. anticoagulated blood. 

A statistical analysis of the variation of the data, comparing fresh 
samples vs. EDTA anticoagulated samples, and for each type of sample the 
variation in viscosity vs. the days of collection, is summarized in Table 
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TABLE 1 
HEMATOCRIT DEP NC I { AND EDTA ANTICOAGULATED BLOOD VISCOSITY 
Standard 95% Conf. . 
Intercept Error Band corr. 


hear Rate = 1.3 sec 


.163 


hear Rate = 0.5 


+42% 


amples 


falling outside the 95% predicted range have been excluded 
from the analysi 


igures in parentheses are minimum and maximum limits 


Paired-t tests were carried out, and 


the significance levels for each of the 


comparisons are given. The analysis was carried out to examine the following 


hypotheses: that there were no differences between the viscosities of either 
fresh or EDTA blood with respect to test day, that there was no variation of 
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TABLI 
EDTA ANTICOAGULATED BLOOD VISCOSITY 


Norm: Significance 
(95%) Level (%) 
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Secon 


t 


rirst 
Second 


Third day 


chelation of the calcium ions by EDTA. There appears to be a fortuitous 
relationship in that at the average hematocrit of 45, the effect of EDTA 
appears to be minimal with much greater effect at lower and higher hemato- 
crits. The importance of recognizing this factor when interpreting the 
viscosity behavior of abnormal blood samples is obvious. Also, the fresh 
blood viscosity shows less variation over time thanthe EDTA samples of the 
same blood, implying that EDTA may interfere with some mechanism tending to 
maintain constant viscosity. 

Limited data on patients suffering from various cerebrovascular diseases 
indicate that the viscosity values are significantly elevated relative to 
what would be expected of the normal range at the given hematocrit. Of the 
16 samples tested, 4 were outside the normal 95% confidence limits, whereas 
on statistical grounds only one would be expected to be. Furthermore, they 
were all elevated. Also, almost all the samples that were within the normal 
95% band were above the best line, indicating that the results were skewed 
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towards higher values. Interestingly, three of the four sickle cell samples 
were within the normal band, although the samples had abnormally low hemato- 
crits. These results would confirm data of previous investigators (10, 11) 
that there are significant viscosity changes in the case of vascular disease. 
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The diarthrodial joints in the human body consist of bones and specialized soft connec 
tissues, each performing complex functions. For example, the articular cartilage in s 


vial joints must sustain high loads on the order of magnitude of several times body weight, 


‘ 


and yet it must provide the very low coefficient of friction during joint motion. The 


ligaments and tendons around these joints must support the joint's stability and 


t 


rapid and high loads. Yet, these soft tissues must not limit joint 1 in functi 


yn 


mn 


oO 


2 


axes. Naturally, the understanding of the mechanical function of these }j« omponents, 


their homeostatic responses, their physiology, and ability to repair re necessary pr 


requisites for the knowledge of joint function. 


The purpose of this symposium is to demonstrate the importance and the for multi- 
disciplinary approach in the understanding of joint function. Fung (1 lescribed tt 
major features of the complex rheological behavior of soft tissues such lung, heart 


e 


e 


, 


, 


urethra, and blood vessels, etc. by showing all these living tissues composite materi- 


als. Different components of the tissue will take up loads at different stress levels 


The loading and unloading of soft tissues follow different and nonlinear paths. The tis- 


sues are anisotropic and possess time dependent viscoelastic properties. urthermore, 
soft tissues have no natural state, i.e., the state a material remembers 1 always re 
turns to when all external loads are removed, to serve as a unique rei for stres 
and strain measures. The soft tissues around the joints also possess e aforemen- 
tioned special rheological features. In addition, effects of aging, physiological wea 
and tear, and the complex repair process of these biological material 

concern to clinicians and surgeons. Articular artilage, f< examp le 

aneural and has little or no ability to repair itself, and yet it must survive for 
time. The healing of a lacerated digital flexor tendon not only volves the strength 
repair, but also its gliding function with peritendinous sheath. 





FIG. 
Mechanical properties of various soft tissues. This diagram demonstrates the 
fiber orientation on the tensile properties. Note that once the fibers are 


slope of the stress-strain curves becomes similar for all tissues. 
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STUDIES OF SOFT TISSUES 


to immobilization and exercise can be best under- 
are known. Akeson, et al (14) demonstrated that 
for a 9-week period resulted in an average 
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ss-link involves an aldol ndensation product that gives rise to an 
unsaturated aldehyde. This is a spontaneous, nonenzymatic reaction 
seems to occur more lily when the collagen molecules are packed 


he helical region of some collagens, parti 
in extensive cross-linking, such as in skin 


Intermolecular cro link form between peptide-bound aldehydes and 
dified-NH» groups of other lysine and OH-lysine residues. This reac- 

is of the Schiff base t . Further stabilization involving reduction 
r addition to the double bond ems to be necessary to render these cross 
links heat- and acid-stable. hes intermolecular cross-links are physi- 
ylogically the most important or 
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FIG. 


Extracellular event leading to fiber formation. 8. Bio 
synthesis of aldehyde near the N-terminal region. 9. Aldol 
condensation leading formation of intramolecular cross 
links. 10. Appearance of aldehydes in the helical region 
of collagen. 11. Formation of intermolecular cross-links 

ff bases). +¢ taggering of collagen molecules with 
in a fiber to generate the 680-A periodicity (640-A in the 
dry state). Arrows indicate the potential sites for inter 
molecules cross-linking (for rat skin collagen). 


ocni 


12. The collagen fibers found in the extracellular matrix of connective 
tissue vary in diameter and distribution. These differences in appearance 
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is very rapid and may occur before the collagen molecule sembles into its 
final position in the fiber. in the other hand, additional ; -+hydes form 
in the helical portion of the molecule during maturation and the average 
number and approximate location of such aldehydes in rat skin collagen ha 
been reported (48). The more insoluble, yet extractable, f 
collagen also contain a significantly large amount of al 
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FIG. 13 


Structures of the most abundant crosslinks isolated from collagens 


after reduction with NaBH,. The asterisk indicates the ite of 


insertion of tritium when radioactive isotope is used for reduction. 
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FIG. 15 


Crosslinks found in various proportions in collagens from 
different tissues. The individual amino acids which give 
rise to these crosslinks, mostly lysine and lysine derived 
aldehydes and histidine, are identified enclosed within 
dotted areas. 


THIOLISM: A DEFECT INDUCED BY PENICILLAMINI 


Administration of penicillamine to animals and humans causes an 
mulation of neutral salt luble collagen in various tissues (51-54). 
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The unusual features of penicillamine posed an interesting question as to 
how crosslinking of collagen could be impaired. As will become apparent 
later, two of the more characteristic properties of penicillamine, namely 
the ability to trap carbonyl compounds and to chelate heavy metals, are of 
primary significance in this connection. The former property manifests it 
self over all the effective dose ranges studied, whereas the latter seems 
to occur only at high dosages, far beyond that administered to humans. 

The mechanisms by which penicillamine inhibits crosslinking of newly formed 
collagen and by which it depolymerizes incompletely crosslinked insoluble 
collagen will be briefly discussed. 

The collagen extracted from tissues of animals treated with D-penicil 
lamine is able to form stable fibers in vitro and is not deficient in 
aldehydes like that from gAP-treated animals. As a matter of fact, its 
aldehyde content is even higher than normal. This suggested that the me 
chanism of action of 8APN and penicillamine had to be different. 

Using the method of equilibrium dialysis, it was observed that cys 
teine, penicillamine and other analogues bind specifically to free alde 
hydes on collagen. Both the free amino and sulphydryl groups are 
necessary for binding to occur. The collagen-mercapto-ethylamine product 
is in equilibrium with its constituents and can be completely dissociated 
by exhaustive dialysis. 


Although all collagen preparations are able to form native collagen 
fibers with native 640A periodicity, only those containing significant 
amounts of aldehydes form stable crosslinks. The collagen from the animals 
treated with high levels of penicillamine is incapable of forming stable 
crosslinks when reconstituted in vitro, and in this connection behaves 
like lathyritic collagen. 


One can conclude from these experiments that penicillamine has the 
ability both to block aldehydes and to inhibit lysyl oxidase activity. At 
lower dosage, it acts primarily by blocking the aldehyde residues present 


on the collagen molecule, whereas at the higher dose it also affects the 
activity of lysyl oxidase (55). This latter observation agrees with the 
fact that 10° * M penicillamine partially blocked lysyl oxidase activity 

in vitro while 10°“ M inhibited it completely . Although the mechanism 
Operating under these circumstances has not been defined, it may be related 
to the ability of penicillamine to chelate divalent cations or to interact 
with other cofactors (Fig. 16). 
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FIG. 16 
Schematic drawing summarizing the mode of action of penicillamine. 
This compound in vivo as well as in vitro will interact with aldehydes 
on collagen, preventing their subsequent participation in the forma- 
tion of intramolecular and intermolecular crosslinks. The reversi 
bility of the collagen defect seen when penicillamine therapy is 
discontinued can be explained by the instability of the thia 
zOolidine complex formed between penicillamine and the peptide-bound 
aldehydes. 
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Thermal shrinkage pattern of strips of dorsal skin (5x50 mm) of 
normal animals and animals treated for 45 days with either 
8-aminopropionitrile (BAPN) or D-penicillamine (0.58 mmole per 
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INTRODUCTION 

nnective tiss f t in the structure of most organs, and the main protein 
component of such C 
which extrude amino ac into a mucopolysaccharide matrix. The basi init formed ir 


> 4} y 


)llagen is produced extracellularly by fibrocytes 
way 1S the three-stranded tropocollagen molecule 


After the formation of tropocollagen, this protein 2S ; in metaboli: 
is inert. The amino acids do not exchange and the macromo > “e ) broken dowr and 
renewed. Consequently, collagen has a property jue 
of collagen has been an ve area of research for many 


T 77 
| 


his work wil 


with only rat tail tendon (hereafter } 
convenient connective tissue specimen from the standpoint of availability id ease 


+ 
sample collection. Experiments with RTT 
properties as early a 954; it has been used in suraery and medi ; before 190f 


were done to ’ acterize ire ind 


bject ¢ t > par 5 te scr e currently ac 
RTT an "a is : tur } nechani< 


Several chemical ibs tances ne se tendon; na re ’ tein cr 
the polysaccharides chondroitin sulfates ; 
cellular materials, and water. The complexity o issue i j increased by a 
hierarchical organization of structures. Thi section describes e structural aspect 

ais t / are presently understood. 


proteoagl yc 


fe) f 


and 
RTT 


Water comprises about 65 if the na > wet weiaq 2 (1.2 jepending on the 
age of the specimen and method of jetermination. Water associates with the collagen 
helix, polar side chains, and the matrix material and is described res; 
structural, bound, or free (3). Four states of water identifiab 


e [ 


freezeable water, freezeable water with heat and temperature of fusion different 


+ fre 


11tTTeren y 


1f bulk water, freezeable water with or the tem 
bulk water, and freezeable water which indistinquishat lk water 
work has been done to understand the ro f water . 
The most important structural omponent o 0 ] the proteir 
represents nearly 99 f the tendon dry weiah (7). Three 104 left- 
chains, two al type and one a2 type, combir ) form a right-handed 
tropocollagen molecule. Tropocollagen has n c ar weight about 
80 nm, and diameter about 1.5 nm (8) 


The primary structure of the 1,952 residues in the al cha j ®) half 
: 
J 


residues in the a2 chain has been solved (9). Little variation equence ¢ 
among tendons from different species, and the salient features p Follows. 
region of the chains yntains 1011 residues, and 9-25 residues terminus are 
helical 


were represented as Gly-X-Y, then proline or hydroxyproline fil pe X¥ perferer 


regions. Every third residue, without exception, is . If the sequen 


The complete residue mposition approximates 33 }lycine, 22% prol yy hydroxyvoroli 
11% alanine, 22% polar residues, and 12% miscellaneous (6). 
hydrophilic, and covalent interactions are possible with such 

7 


these interactions are partly responsible for the hierarchical 
molecules. 





RUCTURES AND PROPERTIES 


he matrix, or ground substi , of RTT is responsible for less than 1% of the 
weight (10). Consequently, characterization of the ground substance lags behind that of 
the collagen. It is known that chondroitin sulfates B and C, hyaluronic acid, and 
proteoglycans are present in the matrix (17). Achilles tendon contains albumin, 
transferrin, haptoglobin, hemopexin, y-globulin, and a glycoprotein called mucoid 2 (11 
Such an elaborate analysis of the ground substance of RTT is not available, but is 
xpected to be similar 
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Hierarchical Oraanization of Tendon (30) 


These tendon components assemble in a multi-level structural hierarchy (12) (fig.1 
The study of structure-property relationships is consequently more difficult and 
interesting. Although the term "collagen" has been applied to each level of organizatic 


+ 


there are structure and proeprty var ions from one level to the next. 


Structure on the molecular level has already been discussed. Essential features 
this lowest level of organization are the single collagen chain with sequence (Gly-X-Y- 
and the tropocollacen triple helix w length 280 nm and diameter 1.5 ma. To form the 
next higher level of organization, fi tropocollagen molecules are intertwined in a s] 
helical twist and are staggered by roughly one fourth their length along the fiber axis 

13). This so-called “microfibri] 1 diameter of 35A and has been reported in 
longitudinal sections, stained cross-sections, and freeze-fractured specimens of RTT (14 


A level only recently included in this scheme is the “subfibril” of 10 to 20 m 
diameter. Eveidence for this classification has come from chemically or mechanically 
disrupted material and from embryonic tissue (14). X-ray diffraction evidence of a 
subfibrillar level has been reported in the case of native RTT (12). 


The next higher entity, the ibril", encanpasses a range of diameters fron 
50-500 nm and is the level of the familiar 65 nm (67 nm wet) periodicity in EM or X-ray 
(14) This banding pattern is thought to be caused by the quarter stagger arrangement of 
tropocollagen molecules within the microfibril and be preserved by appropriate packing at 
the higher structural levels. Fibrils also show crimped waveform. The crimped fibrils, 
together with the ground substance, form the next larger assembly, the "“fascicle”. 
Fascicle diameters between 80 and 320 um are common with a distribution maximum for mature 
RTT at 200 to 240 » (14). There is approximate cylindrical symmetry in the packing of 
cirmped fibrils into fascicles. Occasional evidence for a reticular membrane around 
fascicles has been found. 


The tendon itself is comprised of one to ten fascicles, though it is uncommon for 
this number to exceed four. Tendon diameters are between 200 and 600 um with a 
distribution maximum about 350 nm for mature RTT (14). Tendon samples are easily removed 
from the rat tail and sometimes are surrounded by a reticular membrane (14,1). 
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Crimping 


Much of the recent work in our laboratory on RTT has involved the study of the 
fibrillar crimped waveform. A survey of collagen from a wide variety of sources has shown 
that the waveform is present in many tissues, but only in tissues which are subjected to 
tensile stress (15). Furthermore, data from muscle contraction studies suggest that only 
the strain levelscorresponding to crimped fibril straightening are biologically relevant 
™m Vivo. 


Despite the ubiquity of crimping in collagenous tissue, the causes thereof are 
unknown. An early explanation relied on the contraction of elastin fibers within the 
collagen structure to generate the crimps (14). However, it is accepted that RTT is 
without elastin. Dale et. al. (7) have proposed that contraction of the mucopolysaccharide 
matrix, in which the fibrils are embedded, bucklesthe fibrils. Their model system, high- 
modulus fibers in a soft-elastic matrix, closely paralleled tendon morphology and 
properties. The major drawback of their model is that it incorporated a much larger 
fraction of matrix than is present in tendon. A “wedge dislocation" mechanism |.as been 
recently proposed, but lack of hard evidence relegates this proposal to only speculation 
(14). 


From Gathercole and Keller (15) we have several features of crimping in biological 
tissues. Crimping occurs only in collagenous tissues which undergo biological strain, and 
seems to be a common phenomenon in biological fibrillar aggregates which include only a 
small fraction of matrix. Crimping is independent of collagen chain type, as the 
(al)o(all)) type tropocollagen is found both in bone, which has no waveform morphology, 
and RTT and skin, which show crimping. In addition, fibril crimping has been observed in 
newborn RTT, so the origins of crimping are somewhere within the foetal development (16). 


{IN 





FIG. 


Optical Microaraph of a Fascicle (120x) with Crossed Polarizers as Shown (17) 


Polarizing optical microscopy has been the most important technique for studying 
the crimped morphology of intact, mature, wet RTT specimens (17) (Fig. 2). Diamant and 
coworkers deduced, from the behavior of the banding patterns with rotations about the 
microscope axis and tendon axis, that the arrangement of crimps is planar rather than 
helical, and that the crimp shape approximates a zig-zag rather than a sine curve. For 
specimens from mature rats, the wavelength of crimping is between 20 and 100 .m and the 
angle between the fibril and the tendon axis is between 15° and 20°. Straining the 
tendon decreases the crimp angle until, at some critical strain, the banding pattern 
diseppears (Fig. 3). Teased and thereby disrupted tendons reveal fibrils with a wavy 
structure corresponding to the banding periodicity in the intact sample (Fig. 4). 


Transmission electron microscopy has been used to verify the conslusions listed 
above (18). By sectioning samples at various angles, the planarity of crimping was 
confirmed. Additionally, the existence of sharp crimp apices was directly observed. 
Scanning electron microscopy of intact, freeze-dried tendons showed little information; 
however, teased specimens were seen as collections. of flattened ribbons of parallel arrays 
of fibrils which had a waveform of the requisite dimensions (19). 


A specially designed experiment allowed additional verification of these ideas by 
x-ray diffraction (20). Splitting of the reflections provided a "micro x-ray topography’ 
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Optical Micrograp! f Tendon Fixed without Applied Tension 
(left) and with 1% Extension to Straighten Fibrils (right) (30) 


Separated fr Fascl > Shows Planar 


technique for the 10 m range. Crimp angles can be determined from the angle of the 
azimuthal split for the situation in which the beam is perpendicular to the crimping 
plane (Fig. 5). Previous workers in this area had chosen to avoid this splitting, which 
they considered problematic, by recording diffraction fram stretched fibers (21). 


Age Effects 

An understanding of the aging process is an important motivation for much of the 
connective tissue research. A problem inherent in such research is to define a 
biological age of a specimen. Biological age has been assigned on the basis of tendon 
swelling (22), fatigue (23), solubility (22), chemical contraction (22), and thermal 
denaturation (24,25,26) experiments. tudies of aging have unfortunately clarified only 
the effects of aging rather than the causes, and without direct application of the results 


to living systems. 


The difference betweeen aging and maturation ought to be clarified. Maturation 
refers to development from birth to sexual maturity, which occurs at 3-4 months in rats. 
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FIG. 5 


Low Angle X-ray Diffraction of Rat Tail Tendon Unit 
Photographed with X-ray Beam Perpendicular to the Crimp 
Plane (left) and an Enlarged Tracing Thereof (right) (20) 


Aging refers to changes in the animal after maturity has been reached. Age-growth 
related data will be discussed with this distincticn in mind. 


Cross-linking of the collagen chains, both intra- and intermolecular, is a well 
accepted effect of aging of RTT (22,26,27). The nature of these cross-links is not 
completely understood. From isometric melting experiments, evidence for aldimine cross- 
links between lysine ahd hydroxylysine residues that become more thermally stable with 
age has been proposed. Chemical reduction reactions have pointed to N-terminal 
aldehydes as cross-link percursors (28). In both of these series of experiments, the 
cross-link density was constant with age but the stability of cross-links increased. 


Much is known about the changes in the hierarchical structure with aging. The 
cellular content of tendon is about 45% at birth, declines to 30% by 0.5 months, 10% by 
1.5 months, and stabilizes about 3% at maturity (2). The 35 nm microfibrils have been 
described as loosely packed in new-born tendon and tightly packed in tendons older than 
0.5 months. Several descriptions of the changes in the fibril diameter distribution with 
aging are available (2,29), one of which is presented in figure 6. After maturity, the 
fibril diameter does not change significantly; however, fusion of two fibrils into one 
irregularly shaped structure has been observed (2). 


The crimp parameters, wavelength and angle, change with maturation and aging. 
Foetal RTT fibrils have crimp angles of 35 to 45° and wavelengths of 12 um (16). 
Progressively smaller angles and longer wavelengths with maturation and aging are 
summarized in table 1. The increase in wavelength corresponds to the increase in tail 
length so that the number of crimps is constant (2000 + 300) throughout the lifetime of 
the animal (17). 


TABLE I 


Values of Crimp Parameters, Length and Angle, for Tendons of Various Age (17) 
Age of Animal Number of Samples 


Age of Animal Number of Samples Wavelength (um) Crimp Angle | 
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STRUCTURE-PROPERTY RELATIONSHIPS 


Because tendon serves as the link between muscle and skeleton, it must have a 
reasonably high modulus and toughness while simultaneously being able to flex at joints, 
absorb shock, and retract when muscles are relaxed. The multicomposite structure of 
tendon suits well these requirements. The mechanical and physical properties of RTT are 
important to the study of implant devices, theories of aging, and understanding of 
connective tissue generally. 


The most frequently reported mechanical test has been tensile fracture, from which 
a stress vs. strain curve is constructed. A camparison of the work of various 
contributors is difficult, because of differences in age, strain rate, the definition of 
zero strain, and the ecuinmentused. Results of tensile fracture experiments from 
Kastelic et. al. (30) are shown in Figure 7. Four regions are obvious on the curves for 
mature RTT; a “toe” region, linear region, and two yield regions. The strain rate was 
8%/minute and zero strain was defined as the point at which the specimen was first pulled 
taut. Diamant (17) has published similar results for which the strain rate was 8%/minute 
and zero strain was defined as the tightest sample position that maintains the native 
banding structure under crossed polaroids. Again four regions can be identified for 
mature rats. The results of Nathan (31) and Morein et. al. (32) are from experiments at 
160%/minute strain rate with zero strain being the point of the first needle deflection 
on the force recorder. Despite the huge difference in strain rate, the values for stress 
and strain at fracture fall within the same range (4.0 to 6.3 kg/mm* at 12 to 14% strain). 
These latter two curves display only one yield region however. Van Brocklin and Ellis 
(33) have reported increases in the linear modulus with increased strain rate. 


Chvapil and Hruga (34) reported the fracture stress and strain of RTT at different 
loading rates. They found that the maximum fracture stress occurred over the loading 
range 600-1200 g/min and the maximum elongation occurred over the loading range 180-600 
g/min. These values of load are not directly comparable to strain rate data, but a 
consideration of the average diameter of RTT and shape of the stress-strain curve shows 
that maximum strength of the tendon occurred between 8 and 16% per minute strain rate 
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Stress-strain Behavior as a Function of Age (30) 


and maximum elongation was achieved with a 2 to 8% per minute strain rate. According to 
their data, a strain rate of 160%/minute would produce values of fracture stress and 
strain for lower than one half these maximum values, which is in conflict with the results 
from Kastelic, Diamant, Nature, and Morein listed above. 


The significance of the "toe" region in the deformation experiment is recognized. 
Up to some critical strain, at which all fibrils are straight, there is a gradual decrease 
in the crimp angle. The linear region of the curves is thought to correspond to the 
elastic deformation of the fibrils. Within both yield regions the deformation is 
irreversible, which suggests that structural damage occurs at these strains. Mechanisms 
involving slippage of microfibrils and loss of cohesion between microfibrils are 
implicated to explain behavior in both yield regions (30). 


Attempts have been made to describe the particular contributions of each structural 
component to the mechanical behavior. Dale (7) treated RTT with hyaluronidase and with 
trypsin; hyaluronidase-digestible material included hyaluronic acid and chondroitin 
sulfates C and possibly B; trypsin digestible material was non-collagen protein. He 
concluded that hyaluronidase-digestable material maintains cohesion within the 200 nm 
fibril and is critical for tendon strength and elasticity, and that trypsin-digestible 
material maintains lateral cohesion between fibrils and is not crucial to the strength or 
elasticity. Earlier work with the same enzymes was not conclusive because of difficulties 
with the precision of stress-strain data(10). 


Water is very important to the mechanical properties of RTT. Recent experiments 
in our laboratory have shown that tendons, when dehydrated to 15% water, can reach 
fracture stress and strain values of over 40 kg/mm* and 20% respectively. The values 
exceed the maximum estimate for individual collagen fibrils of 15-30 kg/mm2 (35). 
Interesting curve shapes and the cdisavearance of the "toe" region are other features of 
the work with dehydrated and only partially dehydrated RTT (fig. 8). 


The linear modulus increases, the area of the hysteresis loop decreases, and the 
maximum stress in each cycle decreases for cyclic loading into the “toe" regions (2% 
Strain) of RTT (Fig. 9). After ten cycles, for example, the modulus has increased by 20 
to 40% (2). Rigby (36) has reported similar results toaether with additional evidence 
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Stress-strain Behavior of RTT as a Function of Water Content 
of the Tissue, from 3 Month Old Rats (unpublished results) 


of sharpening of the x-ray diffraction pattern, and explained these results in terms of 
improved orientation. Periods of rest aftercycling, for as long as fourteen hours, showed 
that the values of maximum stress and the hysteresis-loop area returned to their initial 


values (2). 


When tendon is deformed into the linear regions, cyclic measurements indicate some 
changes in mechanical properties, but these are reversible after periods of rest. Crimp 
parameters do not change with cycling. Significant changes resulting from damage to the 
tendon structure do result from cycling into the first yield region, the most important 
of which is the decrease in linear modulus. Cyclic deformation into the first yield 
region also causes re-crimping, which is the appearance of smaller crimped waveforms 


superimposed on the initial waveform. 


Rigby et.al. (37) conducted stress-relaxation measurements with RTT from 4-5 month 
old rats and plotted the data in terms of percentage change of stress as a function of log 
time. For RTT deformed and held at 3.5% strain, the stress leveled at 15% of the original 
value after 24 hours, with a linear region between 0 and 60 minutes and a second linear 
region after 60 minutes. RTT deformed and held at 7.5% relaxed to zero stress after 24 
hours and the curve shape was different. Temperature changes affected the 3.5% samples in 
a reversible manner, but any change in stress level with changing temperature was permanent 
for the 7.5% samples. This relaxation behavior was not variable in the temperature range 0 


OS WG. 


Hooley and Cohen (38) have modelled creep deformation for human digital flexor 
tendon with the assumption that the low-strain viscoelastic behavior of tendon results 
first from shearing of the matrix and then stretching of the straightened fibrils. By 
invoking a crimp angle distribution, they were able to model their experimental data. 
general features of their results are that plots of isochronal creep strain and creep rate 


The 


against stress bend downwards as stress increases. They observed that the non-linear 
viscoelastic behavior of tendon is opposite in several ways to that found for conventional 
polymers. We know of no such study for RTT, but would anticipate very similar results. 
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FIG. 9 


Typical Cyclic Stress-strain Hysteresis Loops at Constant 
Maximum Strain for Deformation into the Linear Regions (above) 
and Yield Region (below) (2) 


Age Effects 


Some change in the mechanical properties is the result of the aging of tendon in 
vivo. The previously mentioned structural changes with age cause mechanical property 
changes, some of which we describe here. In figure 7 the increase of apparent modulus in 
both yield regions with age is depicted. These changes are due to increased inter- 
fibrillar interactions via-cross-linkina, which hinder the principal damage mechanism, 
fibril dissociation. Additionally, it is important to notice that immature RTT shows only 
one yield region, and that the modulus of the linear region remains constant with age (30). 
These data are in accord with previous work (17). The effect of aging on the modulus is 
reasonable since deformation in the linear region is dominated by the collagen fibrils. 
Changes in fibril diameter which occur during maturation are responsible for the increase 
in modulus; after maturity, no fibril diameter increase occurs and consequently the modulus 
does not change. 


Aged tendon responds to enzyme treatment in the same way as 3 month old RTT, but to 
a lesser degree (7). There is a loss of polysaccharides, especially hyaluronic acid, with 
aging and the remaining non-collagen components become more resistant to extraction. 


Modelling the Stress-Strain Nonlinearity 





Modelling of the stress-strain behavior has been almost entirely concerned with the 
"toe" region and crimped waveform. Early attempts have come from workers in the textile 
industry who concentrated on the behavior of single crimped fibers. Other models have made 
use of various non-morphological structural assemblies. Viidik (39) modelled tendon 
behavior with springs, dashpots, and slack elements. Lovin (40) reasoned that the waveform 
is induced by contraction of elastin fibers within the tendon; RTT is devoid of elastin 
fibers. 


Diamant et. al. (8) defined a planar zig-zag with infinitely rigid apices so that 
the concepts of contilever bending could apply. Their results were quite good except in 
the case of high loads or young animals. Cominou and Yannas (41) considered both single 
isolated sinusoidal fibers and a canposite of these. The fit was fair for isolated fibers 
but worsened when the composite was considered. 
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The Mechanism of Resistance to Extension in the Sequential Straightening 


ading Model (14) 


Kastelic (14) has proposed two models which assume a blunted zig-zag waveform. The 
first attributes the stress-strain non-linearity to elastic shear in the matrix. The 
second and more attractive proposal is that a range of crimp angles is present in tendon, 
and therefore, greater numbers of fibrils became straightened and bear load as the 
experiment progresses (figure 10). rimped fibrils are assumed to have no resistance to 
strain. There is some evidence for a crimp angle distribution, and the closeness of the 
model fit makes this a reasonable proposal. 


In the interest of the concise presentation, many details about RTT have been 
omitted. Crimping has been treated preferentially herein, as it has in recent literature, 
because of its significance to the mechanics of the living systems. The thrust of current 
research in this area is in applying in vitro data to the living system. One key to this 
comparison is the influence of water as a transport agent and property modifier. Another 
key is to further characterize the non-collagen components of tendon for a better 
understanding of the enzymatic processes which occur. Finally, emphasis is shifting to the 
causes of aging rather than the effects. 
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ABSTRACT The purpose of this paper is to review the present knowledge of the effects of 
immobility on fibrous connective tissue and to outline present theories about 
the pathogenesis of joint contracture. A severe disordering of connective tissue 
homeostasis appears when the diarthrodial joints are subjected to prolonged im- 
mobilization. Not only are bone, muscle and cartilage affected by disuse, but 


the fibrous connective tissue structures are functionally disturbed. We found 
that the tissues lose significant lubricating and buffering volume of water and 
glycosaminoglycan in concert with increased inter- and intra-molecular cross- 

links of collagen which resulted in joint stiffness. The potential appli 
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FIG. 8 
Radioactive elution profile of a 3N p-toluene sulfonic acid hydrolysate of [~H] NaBH 
reduced collagen from control and immobilized periarticular connective tissue from rabbit 
knee. The abbreviatiations for the radioactive compounds are DHNL, dihydroxynorleu 
HNL, hydroxynorleucine; N&-hexlys, N€-hexosyllysine; Al-His, Aldol-histidine; DHLNL, 
hydroxylysinonorleucine; HMD, Hydroxymerodesmosine; HLNL, hydroxylysinono sucine; HHMD, 


histidinohydroxymerodesmosine. 


FIG. 9 
Elution profiles of the aldol-histidine - dihydroxylysinonorleucine peak from the control 


and experimental periarticular connective tissue collagen rechromatographed on an extended 


basic column. 
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For sof port ermeable biphasic ma ials iiffus f permeation 
ficient t substan al compactior f h é , Fig.4. Thi ompact 
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ABSTRACT Although the internal structure of cilia on micro-organisms and 
in the lung is identical the mechanical principles of fluid 
(mucous) transport may be vastly different. Cilia on micro- 
organisms are generally quite long, less densely packed on the 
epithelium, have a pronounced aperiodic beating pattern and a 
highly coordinated metachronal wave compared to the short 
densely packed, multicellular coordination of cilia found in 
the lung. A simple hydrodynamic theory using two layers of 
Newtonian fluid of differing viscosity is developed to model 
muco-ciliary transport. Features include estimates for the mean 
force distribution, the likelihood that a cilium is located at 
a given height above the epithelium for an idealised beat 
pattern, and allowance for the retardation of the serous layer 
by the cilia. Initial calculations suggest that if the cilia 
just penetrate the upper, much more viscous, layer that mucous 
transport rates are substantially enhanced. Other features of 
the results are the significant retardation of flow in the serous 
layer and the ‘plug’ flow in the mucous layer. 


INTRODUCTION 





The internal structure of cilia on micro-organisms and in the lung are 
identical (Sleigh (18)) showing the well-known 9+2 fibril complex associated 
with cilia and flagella. The cilia bend due to the relative motion of th 
fibrils (Satir (17)). Although the structural and functional (the movement of 
the micro-organism or fluid) aspects of cilia on micro-organisms and in the lung 
are the same, the actual mechanisms of transport may be vastly fferent. 
tris paper we present several theoretical concepts which support this view. 

Morphologically the cilia lining the airways of the lung are different 
from those on micro-organisms. The cilia in the lung are sf often around 
5um. in length, while those on fast swimming micro-organisms re Sum. long. 
The frequency of beat appears to be similar in both circumstance: ] 
although it has been observed by Travani (9) that the cilia be 
the lung increases in the cranial direction. 


On internal mechanical grounds Sleigh (18) suggests that cilia have 
a minimum radius of curvature around 2um. Thus, the longer cilia on micro- 
organisms can develop a very asymmetric beat pattern characterised by the 
familiar pendular effective stroke and the “limp" recovery stroke near the 
epithelium surface. There are excellent hydro-dynamic grounds why this 
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asymmetric beat is desirable and these have been outlined previously by 

Blake (1), Blake et al. (4,5) and Sleigh (18). On the other hand the shortness 

of lung cilia probably limits the beat pattern to sharp pendular motions in both 
the effective and recovery strokes, something similar to that depicted by 

‘heung et al. (6). In the lung, cilia are densely packed on the epithelia 

with only a few diameters separating adjacent cilia (6-10 cilia/ (um) “) 

compared to the relatively spacious arrangement on the fast swimming Paramecium 

(0.1-0.25 cilia/(yum)*) Machemer (14), Blake et al. (5). 


Coordination of the cilia on micro-organisms is well developed into 
metachronal waves whereas in the lung the coordination is more on a multi- 
cellular level. That is, the phase change may be from cell to cell rather than 
from row of cilia to row of cilia found in micro-organisms. Iravani et al. 

(10) observed patches of metachronal fields in rats although they were very 
irregular in size and shape varying from a few cells to hundreds. 


Mucus is secreted by mucous cells. They are abundant on the epithelia of 
the trachea and larger bronchial airways, but are much less frequent in the 
smaller bronchioles, disappearing altogether in the terminal bronchioles 
(Sleigh (18)). At the terminal bronchioles level other secretory cells exist 
producing a watery liquid similar to that in the serous fluid and in the 
osmiophilic liquid lining the alveoli. 


Iravani et al. (10) have reported that mucous transport initiates in the 
form of small flakes, 4-70um. in diameter. These small flakes then aggregate 
to form plaques and then streams of mucus. These streams possibly coalesce 
further up the lung to form the mucous blanket so often reported in the 
literature (Lucas et al. (13); Sade, et al. (16)). Muco-ciliary clearance 
rates vary widely depending on location and circumstances. 


Recent studies of muco-ciliary transport systems (Winet (19) and Gilboa 
et al. (7)) suggests that mucus increases transport rates and in the case of 
the frog palate is essential for transport. When the mucous supply is depleted 
transport is minimal and can only be regenerated by the addition of more mucus 
9x a substance with polymer-like properties. Studies of the ultra-structure 
»9f the mucous blanket by Reissig et al. (15) support the view that cilia 
penetrate the blanket. 


In this paper we use a simple theoretical hydrodynamic model using two 
layers of differing viscosity. If the cilia just penetrate the upper much 
more viscous layer it will be seen that flow rates can be substantially 
increased. We would suggest that one explanation why the cilia in the lung 
are so short is so they can sustain these heavily increased loads that are 
required in the transport mucus. In normal, unloaded circumstances the 
cilia are probably acting sub-optimally, and essentially just “wafting” around 


in the serous fluid. 


Most theoretical models of ciliary transport have restricted their atten- 
tion to micro-organism locomotion in Newtonian fluids. A review of the various 
models can be found in Blake « al. (4,5). Suffice to say that two models 
have been developed to represent ciliary action; in one, the envelope model 
replaces the cilia by an “enveloping” surface over the tips of the cilia, the 
other attempts to incorporate the whole length of cilia within the sublayer 
and hence is called the cilia sublayer model. The theory for the cilia sub- 
layer model has been further developed by Keller et al. (11) and by Liron et 
al. (12). In their models the action of the cilia has been replaced by either 
a volume or a line distribution of force. The replacement of the ciliary beat 
by a line distribution of force is more physically and geometrically reasonable 
than the volume force representation but this gain has to be balanced against 
the (mathematical) analytical simplicity of the volume force representation. 

An extension of the Keller et al. (11) model to an ‘active porous medium’ 
approach was suggested by Blake (3). This is perhaps more appropriate for the 
case when cilia are densely packed on the epithelium as occurs for example in 
Opalina and in the lung. We will develop a hybrid model in this paper. 
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FLUID-MECHANICAL MODEL 





In the mathematical model, we restrict the geometry to two dimensions 
(x,y). We will suppose we have two Newtonian fluids; a serous layer of 
viscosity uw, and depth h, and on top of this a mucous layer of viscosity u 
and depth (H-h) (see Fig. 1.). Thus, the depth of the combined mucous and 
serous layer will be H. The depth of the cilia sublayer (the layer in which 
the cilia operate) will be taken equal to 2. In the analysis, we will con- 
Sider two cases (i) 2 < h, when the cilia do not penetrate the mucous layer 
and (ii) h < 2 < H when the cilia do penetrate into the layer. Actually from 
the mathematical viewpoint, we may ignore the region 2 Sy SH as this moves 
in a ‘plug’ flow with velocity U(2). 


Reducing the mathematical analysis to the simplest level we use only the 
mean-field equation with a mean volume force f dependent only on the height y 
(see Keller et al. (11)). Thus, the relevant equation in the serous layer 
(0 Sy Sh) is 


(1b) 


"2 dy 


where y, and yu. are the viscosities in the serous and mucous layer respectively. 
The boundary conditions and matching conditions at the interface between the 
two Newtonian fluids are, 


u (0) = 0 
u (#) =O | 
U (h )= U (h) 
My u (h ) = 


where the dash implies differentiation with respect to y. 
Equations (2c) and (2d) imply continuity of velocity and stress across the 
interface respectively. 


In presenting the solution to these sets of equations it is much 
clearer if we separate the cases where the cilia do and do not penetrate 
the mucous layer. The solutions are formally presented in terms of integrals. 
For the case when the cilia do not penetrate the mucous layer (2% < h), 


K(y,&) £ (&)dé (3a) 


U 


where 


In this case the mucous layer moves with the same velocity as 
of the serous layer in a plug flow. 


For the case when the cilia do penetrate the mucous layer we 
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FIG. 1 


(a) Illustrates geometry and arrangement of layers used in the mathematical 
model. (b) Idealised ciliary beat pattern used in analysis. 


In a single Newtonian fluid Wu (20) has suggested that one optimal solu- 
tion for a given rate of working allows for the force f to vary as sin 
We can do likewise in this double fluid model with the additional condition 
that the force distribution is also linearly proportional to the viscosity. 
This implies there will be a discontinuity in the force distribution at the 
interface. If we substitute this sinusoidal force dependence into both 
equations (3) and (4) we find for the case when §8=1000 and the penetration 





depth into the mucous of 0.12 that the flow rates are increased by a factor of 
50 for the case when the cilia penetrate the mucus. This suggests that the 
penetration of the cilia into the mucus in this model is the predominant 
mechanism for transport. This contrasts with the action of cilia in propellin 
micro-organisms where the pronounced aperiodicity and asymmetry of the ciliary 


beat between effective and recovery stress is the prime mechanism for 


This model also suggests a massive increase in the rate of working of 


when they penetrate the mucous layer. Presumably the cilia are capable of 
coping with this increased work load. 


However there is a major inconsistency with this model. 
flow rates we calculate an almost constant shear zone (Couette 
serous layer. This is the main weaknesses of this model as it 
for the additional resistance of the cilia due to the higher 


SUBLAYER RESISTANCE MODEL 


For any realistic model we need to account for the additional resistar 
of the cilia due to increased flow rates attributed to the much more 
upper layer. The force distribution acting on a cilium is directly 
tional to the difference between the cilium's velocity and the local fluid 
velocity. To obtain an average volume force distribution we need to 
over both space and time for a given ciliary beat cycle 
(11)). The aim of this paper is to illustrate, in a qualitative manner, ths 
important physical processes involved with muco-ciliary transport. To do tl 
we will make certain assumptions about the volume force distribution and 
ciliary beat cycle. We will suppose that 


V1 
propo 


averace 


(see e.g. Keller et 


where U (y) is an averaged cilium velocity, “(y) 15s the average over a 
the cosine of the angle between the normal to the cilium and the x-axis 
function of height y and is a cilia density parameter with dimensions 
inverse length. In the next section a method of deriving both 

E(y) for an idealised ciliary will be presented. 


IDEALISED CILIARY BEAT PATTERN 


We will suppose that the effective stroke of the ciliary beat 
consists of a pendular clockwise rotation from to 8 ( 


see Fig. 
cilium length will be taken as equal to £ with 


radian frequency ; 
in the recovery stroke, we will suppose that we have an anti-clockwi 
from 8. to 6: the effective cilium length now being 2 and with 
frequency : 
then, 


radian 


r : 
The cartesian coordinate expressions for the beat pattern 


Effective 


overy stroke 
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Here s is the length along the cilium. It is convenient if we non-dimensionalize 
the linear dimensions with respect to 2, and velocities with respect to ff 
(where f is the frequency of the Ciliary beat). 


This non-dimensionalization 
produces 


Y = y/ U* = ’ U* U /fi 

Q £2 c c 

and T= 22. The force distribution on 

the formula from the appendix of Blake 
the following expression. 


the model cilium is obtained by using 
(2) and by averaging over time we obtain 


obtain this expression we have set L.. = Lsin 8) where L. = he and 
2 /2. We would obtain slightly different profiles if we used either 
: e 
inequality, but the profile is simplified by assuming equality. 


The three layers of the profile correspond to (i) the region in which the 
cilia rotates in both effective and recovery stroke. To derive this expression 
we are considering the time integral of the force and hence the velocity because 
of the linear dependence between the two. However, since we return to the same 
location, this integral must be identically equal to zero. (ii) In this region 
the cilium is always present in the effective stroke, but only for part of the 
recovery stroke. (iii) Only the effective stroke operates in this region. 


We can also obtain €(Y) for this special beat pattern. 


pH , G- aApsin 20), 


oO 


, 2¢ 2/y (l- 1/y)sin26, 
© 2(1+w Jw) © 2(6 =8,) 40 7u 
= oO l oe f£ 





(1+ L/y) cos (¥/Le) , (1- 1/y) 
(6 -6,) (1l+w /w_) (6 -6,) (1l4+w /w 
é r re) l e r 





oO l 


To calculate © we have neglected the mean of the time dependent components of 


the velocity and the orientation tensor. However €(¥) recognises the fact that 


the cilium is only located at a given height for only part of its effective and 
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recovery strokes and hence is very important in determining the resulting 
velocity profile. In Figure 2(a) graphs of U.*(Y) and €(Y) are shown for the 


above idealised beat pattern. 











FIG. 2 


The functions U* and © as defined by (7) and (8) respectively. To 


obtain this plot we have taken w /o. = 3.0. 
e 


NUMERICAL SOLUTION 





To develop a better model we solve, 


T (us (¥) - E(y)u(y)) (9) 


numerically for the values of us and € given in (21) and (22). 


Actually in the region 0 SY S L_sin 6 with U* = 0 and € 
obtain an analytic expression for the velocity 


U = U sinh T € ty (10) 
° ° 


with _ an unknown being determined from the numerical solution at Y=L sin 


If T » 1 we might expect that the ‘active porous medium’ result of Blake 


u(y) us (¥) E(¥) (11) 


will be a good approximation. 





importance of 


© estimate the numerical range 
umed that y, the ratio of normal to 
is two (in reality it will be omewhat 
estimate for from the normal resist 


ylinders (see Happel et al. (8)) as f 


radius of the cilium 
consider examples at 





Active porous 
medium 





solving (23) numerically for (a) 5 and 
1,10,100 and 1000) are indicated on the diagram: 


active porous medium’ velocity prof > given by 


liagrams. 


c 


» and plotted the velocity profile for 
corresponds to the case when the lower and 
sity or when the cilia lie entirely within the 
increasing the velocity of the fluid increases until it 
approaches ost 1e maximum velocity for 6 = 1000. In this example the 
flow ra ‘ the mucous h increased by an order of magnitude for 6 1000 
over f l. Thus penetration of the cilia into the mucous layer substantially 
increases the flow rate in this example. Other important features to notice 
are the retardation of the velocity in the serous layer by the cilia and the 
almost constant ‘plug’ flow the mucous layer. Also in Fig. 3(a) we have 
included the profile corresponding to (11). Fig. 3(b) is similar to Fig. 
3(a), but with two enhanced features of the velocity profile. The first 
concerns the much larger velocity of the mucous layer in the case of 6 l, 
such that increasing the viscosity of the upper layer a thousand fold only 
doubles the velocity. The second feature is the almost total retardation of 





flow in the serous layer, in fact it closely approaches the ‘active porous 
medium’ profile which is also included in this figure. 


This example suggests that the penetration of cilia, 
important for enhanced flow rates, does not produce the pro id increase in 
flow rate. It seems that cilia density on the epithelial surface may play 
important role in the efficiency of mucous transport. 


ar 


CONCLUSIONS 

The theoretical models suggest that the penetration of cilia 

mucous layer enhances transport, but the extent of the enhancement 

the subject of debate. Two other features that are important are 

flow in the mucous layer and the retardation of mean flow in the 

To develop such a simple model many assumptions were used in the analysis 

although all the essential physical features of muco-ciliary transport were 

except for the very important aspect that the fluid in the mucous 
= 


»n-Newtonian. The memory effects of the mucus may well provide 
another important mechanism for transport because of the aperiodicity of 
ciliary beat. To study this analytically, even for linear models, w 

involve a much more detailed analysis of the full unsteady equations - motion. 
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ON THE MECHANICS OF MUCOCILIARY FLOWS 
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ABSTRACT Cilia-generated fluid flow between two strips of frog 
palate epithelium has been measured and flow velocity 
profiles have been obtained from the data. The strips are 
mounted in a covered chamber so as to form a ciliated 
channel and fluid motion is marked by suspended 1.0 um 
polystyrene latex spheres which act as tracers. Results 
from a control non-mucus system and a 40% - 50% v/v mucus 
system are compared with predictions of a theoretical 
model. This theory for flow in mucociliated channels 
incorporates flow due to a static pressure drop between 
the channel ends and a difference in the cilia tip speeds 
between the two strips. 


INTRODUCTION 


Cilia which propel liquids in tubes of mammals are found in the 
otolaryngeal (e.g. Eustachian tube), respiratory (e.g. tracheal tube), 
nervous (e.g. brain ependyma) and reproductive (e.g. oviduct and 
ductuli efferentes) systems. In the Eustachian tube cilia help pre- 
vent large changes of pressure in the middle ear by clearing out 
liquid plugs, (1). In the trachea cilia remove mucus and particulate 
matter. In the ductuli efferentes cilia may keep spermatozoon sus- 
pensions moving toward the epididymis and in the ampulla of the 
oviduct they are the main source of propulsion for the ovum (2). In 
the brain ependyma cilia enhance transmural transport (3) by reducing 
the unstirred layer. 


Morphologically the cilia lining biological tubes such as the 
trachea, cervix and oviduct are different from those of protozoa. 
Typically the tubal cilia are short, usually around 5 um in length 
while those of fast swimming protozoa are 7 - 15 um long. The beat 
frequency appears to be similar in both circumstances, 10 - 30 Hz, 
although it may change with time, as a function of estrus cycle phase, 
(4) (however, see also (5)) and location (as a function of distance 
from the cranial end of the lungs (6)) in tubal epithelium. 
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Cilia of ciliated epithelium are usually crowded into tufts on 
cells which may be separated from other ciliated cells. In the lung 
where ciliated cells form continuous laygrs only a few diameters 
separate adjacent cilia (6 - 10 cilia/um¢). In the oviduct where 
ciliated cells comprise 40 - 80% of surface cells, and carry 200 - 
300 cilia each (7) the spacing is of the order of 0.05 to 0.15um or 
8 cilia/um“ (measured from (7)). In contrast the spacing of cilia 
over the surface of the fast swimming Paramecium is 0.1 - 0.25 cilia/ 

(2 9). 


\V> 


Ihe fluid propulsion generated in tubes is termed “internal flow" 
because the liquid is bounded on at least two opposing sides. There 
are well established fluid mechanical grounds for expecting that 
cilia-generated internal flows are different from “external flows" 


(flow in fluids bounded on one side only) particularly with respect to 
the existence of a pressure gradient along the internal system. 
Accordingly, we cannot normally use external flow systems as experi- 
mental models for internal flow models. 


In mucociliated epithelium of both the respiratory and reproduc- 
tive tracts “secretory cells" which contain glycoprotein granules and 
Clara cells are interspersed with ciliated cells (10) in varying 
degrees. The concentration not the oligosacharide ratio of glycopro- 
tein suspensions secreted from secretory cells appears to be the main 
determinant of its rheology. Thus low concentration (< 0.5% w/v) 
Suspensions with Newtonian rheology are secreted from serous secre- 
tory cells and high concentration (> 0.5% w/v) glycoprotein suspen- 
sions with non-Newtonian rheology (at lower shear rates) from mucous 
cells. The glycoprotein secretions from mucous cells of the cervix 
will reach concentrations as high as 7.5% (w/v) (11), and during 
ovulation low enough such that the fluid can be considered Newtonian. 
In the lung and oviduct the secreted glycoprotein concentrations cover 
a wide range simultaneously. At one extreme they are great enough to 
form a blanket layer (12, 13, 14) over the cilia and at the other 
extreme they form a watery serous layer amongst the cilia. The tran- 
sition to the mucus blanket state apparently begins, in the lungs at 
least, when small aggregates or flakes of mucus gel, 4 - 70 um in 
diameter are transfered from their secretion sites to points of 
aggregation where they form larger plaques and then streams of mucus 
(15). These streams possibly coalesce downstream into blankets. 


The ultimate internal flow experimental model would be the muco- 
ciliated tube. But development of merely ciliated tube biological 
models has been hampered by an inability to resolve the flow paths of 
tracer particles in the lumen. So-called transparent tube walls are 
Simply not transparent enough to allow the special illumination 
optics (e.g. darkfield, Nomarski interference and Hoffman modulation 
contrast) which are needed to resolve small tracer particles. Fiber 
optics in conjunction with “windows" in the tube wall have improved 
resolution considerably; but the fiber interferes with flow. 


In contrast, flow in ciliated channels places no optical obstruc- 
tion other than the cilia (or "“mucocilia") between the observer and 
the light source so fluid tracer particles are easily illuminated. 
Moreover, it has been shown by Liron (16) that ciliated channel flow 
velocity profiles are qualitatively similar to those of ciliated 





tubes. Accordingly, the ciliated channel is a reasonable qualitative 
approximation for ciliated tube flow and a significantly more correct 
experimental model than the external flow system. 


We report here observations of cilia-generated flow in both 
ciliated and mucociliated channels. Both sets of observations are 
compared with predictions of flow velocity profiles from fluid mechan- 
ical theory described below. The experimental models are constructed 
with frog palate epithelium. 


METHODS 


) 


Large (> 7 cm body length) Rana pipiens and R. catesbiana are 
maintained in plastic containers with ca. 1 cm deep distilled H20 at 
16 - 17 °C in an incubator. Frogs are removed as needed and immedi- 
ately processed on a cervical dislocator (Cervical Dislocators, Inc.). 
The dorsal part of the head is then excised and immersed inc. 20 ml 
amphibian phosphate buffer solution APBS (17) until needed. The 
mucus is collected by mounting the severed end of the excised head on 
the open end of a graduated centrifuge tube and allowing it to drain 
for about 10 minutes. 


The epithelium is then prepared by carefully excising it from 
the palate with microsurgical scissors in one or two (first cut along 
the midline) pieces. The pieces are then placed cilia-down on APBS- 
moistened filter paper and allowed to ‘walk off' excess mucus. This 
walk-off also serves as a check for cilia motility. Ciliated strips 
are obtained by gently teasing the epithelium into a flattened 


position and cutting longitudinally through both the epithelium and 
filter paper. Both length and width of strips vary but a minimum for 
each of 5 mm is maintained. Maximum length does not exceed 15 mm. 
During and after cutting the epithelium is kept moist with APBS. The 
freshly-cut strips are then gently drawn onto a coverglass panel to 
which they are tied with nylon (#50) thread. 


Two ciliated panels are then placed in the observation chamber 
as shown in Figure 1. The chamber contains either a control solution 
APBS+, which is APBS with less than one drop/20 cc of | um polysty- 
rene latex spheres added, or APBSm which is APBS+ mixed with the 
drained-off mucus collected earlier. The various dilutions of APBSm 
are obtained by adding APBS+ and mixing vigorously with a wooden rod. 


As indicated in Figure 1 the observation chamber fluid is contin- 
uous with a fluid reservoir and completely surrounds the ciliated 
panels such that the midpoint of the chamber between the two ciliated 
panels is at least Imm from a rigid wall. The ciliated channel width 
which is the distance between the parallel ciliated panels is adjusted 
by simply sliding the panels. The chamber construction also allows 
the alignment of the panels so that their cilia may beat in opposite 
directions. 


Flow of the polystyrene fluid tracers and the various ciliary 
beat parameters are recorded on motion picture film through a Nikon 
microscope (LKe) utilizing Hoffman Modulation Contrast optics and 
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strobed—by camera—xenon arc illumination. Focus during filming is 
maintained in the vertical midplane. Where channel width exceeds the 
microscope field, overlapping field sequences are filmed at each 
ciliated panel and the channel width is measured at a lower magnifi- 
cation with a calibrated Whipple disc. 


Each series of observations consists of filming at various 
channel widths with the widths being varied in a random fashion. 
Random channel width variation is necessitated by a possible relation- 
ship between time following excisement and ciliary activity as 
indicated by the generally short duration of ciliary activity in our 
chamber as compared with the ciliary activity of an uncovered epithe- 
lium still attached to the excised dorsal head (cf 18). 
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Frog ciliated epithelium biological model. 


The films resulting from the above observations are analyzed 
with the aid of a NAC film reader. The fluid tracer motion of the 
polystyrene spheres is plotted frame-by-frame and the average velocity 
of a particle or group of particles at a given distance from the ep- 
ithelium determines each point of a flow velocity profile plot. 
Where ciliary beat frequency f cannot be resolved directly, and wave 
velocity c and wave length \ can, f is calculated from 


Cc f) 


The Fluid Mechanical Theory 


We suppose the cross-section of the ciliated channel is as illus- 
trated in Figure 2. The distance between the tips of the cilia sub- 
layer (the extra-ciliary width) is D'=2b (=D - 22). The depth of the 
chamber is much greater than D' so we can neglect the variation in the 
x direction. The velocities generated by the cilia at the top of the 
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sublayers are U; and U, in a direction out of the page (i.e. in the 
z direction) lihen mutus is present it is assumed to be completely 
penetrated Be the cilia and so is taken to have a layer depth of 
g-5 where 6 is the thickness of the serous fluid. 


In the following analysis all fluids are taken to be Newtonian 
with a known density and viscosity » although we can consider the 
fluid in the layer closest to the epithelium (the serous layer) to 
have a lower viscosity as previously discussed by Blake & Winet (20). 
‘Because of the low flow rates and small dimensions the Reynolds num- 
bers are very small, allowing us to apply the Stokes flow equations 
in one dimension 


where p is the pressure, and U the fluid velocity in the z direction. 
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FIG. 2 


The theoretical fluid mechanical model 


The flow-field can be divided into two parts; (i) the two cilia 
Sublayers and (ii) the central lumen region. Any pressure gradient 
will have an extremely limited effect on the fluid velocities in the 
cilia sublayer because of the increased resistance due to the surface 
area presented by the cilia. Accordingly, we will suppose that in the 
cilia sublayer the flow is due entirely to the beating of the cilia. 
In a previous paper (20) we described a simple theoretical model using 
serous and mucus layers of vastly differing viscosities to model as- 
pects of mucociliary transport. In this report it is sufficient to 
approximate the flow field in the cilia sublayer as follows, 





-(b+2). 


That is, the mucus section (of depth 2-6) of the ciliary layer moves 


along as a blanket with a constant velocity Um (Uj or U2), while the 
serous layer is presumed totally retarded by the dense mat of cilia. 
This is probably a reasonable approximation for channels in chambers 
where chamber depth >>D as is the case with the present model. Such 
a conclusion is suggested by our previous theory (20). (There is 
evidence a damping of fluid flow by the glass walls 

that the approximation is not valid for chambers with thicknesses of 
the order 2(b+2) or less (19).) 


For the non-mucus control we have to match the sublayer flow with 
the external flow in one continuous profile but since our experimental 
model shows only the extra sublayer flow, we shall present only this 
part of the flow equation which is the same in both the control and 
mucus models. 


In modelling flow in the central lumen we divide the velocity 


field into two components induced by (i) the pressure drop along the 
apparatus and (ii) the cilia 


(4) 


2 
D' ld 2 2 
) =a ge ly - be) ~~ (4b) 
The volumetric flow rate per unit length in this system 


0 $03 (8) 


f° 
- )(Uy+Up)+ J Uly)dy = 2(2-6)(U,+U, 


We can usually neglect the first term on the right hand side of 
(5) as it is small in comparison to the other terms; and in the fol- 
lowing analysis it will be dropped. To uniquely determine the flow 
velocity profile either Q or the total pressure drop per unit length 
\P/L = YP (where L is the length of the sealed chamber) must be pre- 
scribed. Once one is known the other can be obtained from (5) by 
using the substitution dp/dz VP. Unfortunately, neither Q nor yP 
was obtained in this experiment but we can use the observed flow 
velocity profile to obtain a measure of YP and at the same time check 
to see if the present level of approximation in the theory is accep- 
table for predicting profiles with a reasonable fit to the data. The 
evaluation of the theory, accordingly, will consist of plotting the 





observed profiles and then using measured Up values to calculate and 
plot U(y) for different values of VP using (4). The VP giving the 
best fit will then be considered the best estimate of the experimental 
model YP. The remaining parameter » will be kept constant at that of 
water although it will be considered as a parameter for the muco- 
ciliary systems where mucus may invade the lumen of the channel. 


Contro1-APBS+ 


Flow velocity profiles for APBS+ in three ciliated channel widths 
D are shown in Figure 3. All values are non-dimensionalized by their 
reference. Here is the thickness of the cilia sublayer and r is the 
extended length of a single cilium. We take =r 5 um for this re- 
port. Also U. is the maximum particle velocity. The measured U, is 
assumed to be the true Um even though the magnifications utilized do 


(a) 


























Control flow velocity profiles for ciliated channels of 
Ordinate is reduced distance from lower epithelium sur 
is reduced velocity. Curved arrows indicate direction 
fective stroke. Arrows at data points indicate flow 
of arrows indicate standard deviation of the sample. 
est U, is used to non-dimensionalize in each case, 


Li 


in text. 





N MUCOCILIARY CHANNELS 


not allow measurement of U within the sublayer. This assumption is 
based upon earlier observations in which Cheung & Winet (19) indicated 
that the true U. lies outside the sublayer. 


In most preparations the ciliary activity on both ciliated panels 
do not match. Moreover, when D exceeds the microscope field and two 
overlapping film sequences are required, it is possible that there is 
a change in ciliary activity during the scene shift. Accordingly, 
even though the profiles illustrated represent the best wet matches, 
one can detect in 3b & c that Up is achieved at one panel only. 


Table 1 presents the various Up's as well as the variations in the 


distance from the panel surface at which particle velocities drop to 
50% of U 
m* 


As one follows the curves toward the channel center from the U 
position a rather rapid decrease in U occurs. The rate of this de- 
crease is greatest in the narrowest channel but not smallest in the 
widest channel as indicated in the last column of Table 1. 


Backflow is present in channels as wide as D/2 = 420 but is ab- 
sent in those as narrow as D/2 = 125. The dashed part of the Figure 
3b curve is an acknowledgement of the lack of tracers at the very 
center of the channel which would have provided more direct evidence 
of the lack of backflow. 


l. 
TABLE 
in um sec 


absent 
18 - 47 


9 - 18 


The point at which flow changes direction (U = 0) is not the 
same on both sides of the D/2 = 420 channel (Fig. 3c). On the side 
with the larger U_, the transition point is y/2 = 80 but on the other 
side (upper part OF the figure) U = 0 at y/2 = 60 (420 - 360). 


Some measurements of U, as a function of f have been obtained 
and are presented along with the assumed point Up, = 0 at f = 0 in fig- 
ure 4. In all cases f is calculated from measured c and A values. 
Um used in each case is the average velocity of the fastest particle 
observed (i.e. not the average for all particles at that distance from 
the ciliated surface; it is felt that the fastest particle is the most 
‘unencumbered’ tracer, other particles may have small strands of mucin 
attached which may be stuck to slower-moving surfaces). 
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FIG. 4 
Flow velocity of fastest measured tracer particles as a function of 


beat frequency f. Plot is extrapolated back to the point (0, 0). U 
ine -6 + cn a m 
is in 10°" ms and f is in Hz. 


Mucociliary Channel-APBSm 





"Vigorous mixing" of APBS+ with drained palate mucus suspension 
does not result in a homogeneous suspension. Strands of mucus are 
still evident in the mucociliated channels particularly when the two 
epithelia are allowed to touch and then separated. Even without the 
aid of Hoffman Modulation Contrast illumination one can see an inter- 
mittent series of mucus strands and patches wafting over the cilia. 
We record no data when the thickness of this mucus ‘layer’ appears to 
be greater than 2. The vigorous mixing does allow a somewhat random 
distribution of mucus patches in the reservoir fluid so that the 
channel has a constant supply. 


Flow velocity profiles for APBSm in three ciliated channel widths 
are shown in Figure 5. The percentage values are in volume of drained 
mucus divided by the total suspension volume. The graph variables are 
the same as in Figure 3 and so are the apparent sources of U_ variation. 
In the narrowest mucociliary channel, however, the difference in mot- 
ility between the two ciliated surfaces is so great than an adjusted 
curve is added to Figure 5a. Of the two curves the solid line is 
based upon the absolute U. (which is Uz) and the dashed line is a 
split of the plot with the upper part based on Uj and the lower part 
based on U The resulting curve which is formed by matching U, and 


U, -based furves is more symmetrical than the plot based on one Ue 


__. The U; and Us values in figures 5b & c match fairly well so a 
splitting of U. is not warranted. 


The various Up's and 50% U-drop distances for the mucociliated 
channel are presented in Table 2 as a parallel to Table |. 
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In general the values of 
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FIG. 


ciliated channels of three width 
two values of U, the solid line is re- 
epithelium ana the broken line to U, 
the true U,). The symbols t, and t, in 


thelia were not photographed Simultane- 


mbols are defined in the text or legend 


TABLE 


sec 


are smaller for ‘mucocilia’ than they are 


wing the profile from U_ toward the cham- 
> general pattern as the Control; i.e. the 
narrowest channel exhibits the greatest slope and the mid-width chan- 
nel the smallest slope. 





CILIARY HANNELS 


One notable difference between the control and mucociliary flow 
profiles is the latter's apparent non-uniform approach to backflow 
with increase in channel width. In the control the minimum particle 
velocity Ug decreases with increase in channel width in both its di- 
mensional form (Ug) and its non-dimensional form Up/Up. In the muco- 
ciliary channel the non-dimensional values (see Figure 5) increase be- 
fore they decrease, but the dimensional values decrease uniformly with 
increase in channel width as shown in Table 3 which compares the two 
systems. 





| APBSm 








FIG. 6 


Reduced flow velocity of the slowest or most negative measured tracer 
particles as a function of channel width. In each case U_/U_ was 
esultings 
values averaged. Error bars are standard error of the sample. Other 


symbols are defined in text. 


calculated for each epithelium of a channel and the two 


CONCLUSIONS FOR EXPERIMENTAL RESULTS 


The differences between ciliary and mucociliary flows in these 
experiments are major only in the values of U, attained in comparable 
channel widths and the drop in particle flow velocity from Um to Uo. 
Since no measurements of f were obtained from the mucociliary systems, 
we cannot be sure that these differences were simply rheological eff- 
ects of the mucus. Moreover, no bulk viscosity measurements were ob- 
tained so the difference in mechanical load upon the cilia of the two 
systems cannot be assessed. We can, however, be reasonably sure from 
the observations of Machemer (8) that the change in load from 40% to 
50% suspensions is not sufficient to cause a significant change in f. 
Accordingly, the APBSm results may be lumped together as a represen- 
tative unit. 





CILIARY CHANNELS 


3. 


TABLE 


D System U/U, 


7.2 0.54 
26 0.46 
124 T 0.62 
125 0.275 
156 -0.30 
420 -0.30 


The generally lower Um values for mucociliary flow are probably 
not the result of a reduction in f because these are the velocities 
at which the mucus layer moves and with the exception of the med- 
width channel compare well (see Table 2) with the 167 um s~' mucus 
layer velocity measured by Gilboa and Silberberg (18). 


+ 
pith ——e 


180- 


D/r = 418 











FIG. 7 


Comparison of theoretical with experimental control models via pre- 
dicted vs. measured flow velocity profiles in ciliated channels. 
Heavy line plots are theoretical. Solid circles are data points from 
Figure 3. Open circles are extrapolated data points. All symbols 
are defined in the legend of Figure 3 er the text. 
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It is evident that at some channel width in the range 124<D/2<420, 
= 0 (the transition point). Beyond this point U,<0. If Uo/U,, is 
a linear function of D/2 then the control plot drawn in Figure 6 is 
valid and the transition point is around D/2 = 270. If U/U,, is a 


nonlinear function of D/2 as suggested by the APBSm plot of the fig- 
ure, then more observations will be necessary before the form of the 
relationship is clear. The data summarized in Figure 6 suggest that 
mucus allows backflow to appear in narrower channels than the control 
does. In the present cases the control channel has to be almost 
twice as wide as the mucociliary channel before backflow appears. 
More data for both regimes will be necessary, however, before any 
conclusion is possible. 
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FIG. 8 
Comparison of theoretical with experimental models via predicted vs 
measured flow velocity profiles in mucociliated channels. Format 


for plots is described in the legend for Figure 7. All symbols are 
defined in the legend for Figure 5 or the text. 


DISCUSSION 


The Agreement Between Theory and Experiment 





The critical test of a theoretical model is the degree to which 
it conforms to the experimental model. Theoretical and experimental 
results are compared in Figures 7 and 8. Here the plots from Figures 
3 and 5 are reproduced with data points and the theoretical results 





from equations (4) are added as solid lines. In some cases dashed 
lines representing better "by-eye" fits are added for comparison. 


In all cases the theoretical curves were generated by requiring 
that they run through three points, Uj, U2 and the minimum value of 
U, Up. In cases where data was lacking for Up the extrapolated 
values from figures 3 and 5 were used. Figure 7 represents a ciliated 
channel propelling a single Newtonian fluid (the control). We might 
expect these profiles to compare best with the data because of min- 
imized error due to the absence of mucus. (Slow-moving strands of 
mucus can pull attached tracers back and/or fast-moving mucus strands 
which have just broken can pull attached tracers forward during their 
rapid recoil.) There is indeed better agreement between the Newtonian 
fluid data and theory than between non-Newtonian fluid data and 


theory which is plotted in Figure 8. The VP's and values of Uo, are 
compared in Table 4. 


7P 
min! Um dynes-cm_ 


Conditions U 


+ 0.515 2.17 x 10 
0.281 3.56 
-0. 398 4.2 
0.478 
0.688 


154 n -0.403 


The variation in data from the mucociliary channel is not limited 
to fluctuatior ibout a th theoretical curve as shown in Fiaqure 8. 
Indeed an auxiliary curve had to be drawn (dashed line) for the D'/2 = 
124 & 156 cases; the latter to allow for the probability that U] & U2 
changed between t] and to. Also, if the point-to-point variation in 
the data characteristic of the mucus channel plots may be attributed 
to mucus strands pulling tracers, the relative flatness of the 
D'/2 = 124 channel profiles cannot. This result has to be due to the 
smaller U, values (see Table 2) which is probably due in turn to re- 
duced ciliary activity. 


One important feature to note in both the theory and experiments 
is the importance of the width of the channel D or D', the viscosity 
of the fluid in the central lumen and the presence of no-slip side 
boundaries on fluid flow. For a given positive pressure gradient 
(dp/dz > 0), the wider ti hannek the weaker the slow an the centrat 
section and in the widest channels of those included in our obser- 
vations the flow actually reverses (see Figures 3c, 5c, 7c, 8c). 





Similarly, the wider the channel the more important the retardation 
effect of the no-slip (slide and cover-glass of chamber) walls will 
be; although in the present model the spaces between the epithelium 
and the chamber walls reduces this factor considerably. However, in- 
creasing the fluid viscosity has an opposite effect, although it may 
well reduce the fluid velocities near the cilia tips and overall in 
the channel as indicated in the comparative values of Tables | and 2. 
An increase in viscosity will reduce the influence of the pressure 
gradient which must now balance the increased viscous stresses. One 
consequence of this prediction is less parabolic central flow velocity 
profiles which is indeed observed as indicated by a comparison of 
Figures 3ab with 5a (split) and b. 
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ABSTRACT A series of experiments have been carried out to investigate the 
muco-ciliary transport in the trachea of rat, rabbit and chicken 
with cinemicrographic analyses of the movement, beat pattern 
and metachrony of the tracheal cilia. As a principal test case 


of the cilia-tip-penetration-into-mucus hypothesis, detailed 
studies have been made on the transport of a biochemically via- 
ble mucous plug artificially introduced into the chicken trachea 
4 
4 


both with and without a pressure differential across the mucous 
plug. Guided by the result of these pilot experiments, a theo- 
retical two-layer model of the muco-ciliary transport is intro- 
duced, in which the cilia are assumed to penetrate shallowly 
into the overlaying mucous layer during the effective stroke and 
to withdraw from the mucous layer during the recovery stroke. 
Both the mucus and the serous fluid are assumed to be 
Newtonian and with widely differing viscosities. Based on this 
model the mucus transport rate is found to depend linearly on 
the ciliary beat frequency and also on the time of ciliary tip 
penetration in the mucus. Results are also given on the pro- 
pulsive force contributed by each individual cilium and an esti- 
mate of shear and shear rates within the mucous layer. 


IDUCTION 


Muco-ciliary systems occur in v us organs in mammals (and m vertebrates) 
such as in the respiratory and reproductive tracts and in olfactory, auditory and other 
sense organs, Although they differ in physiological functions and there is considerable 
morphological variation in different organs, the muco-ciliary systems in the respiratory 

. 


and reproductive tracts are essentially for the purpose of transporting biological and 


extraneous matter. The present study is devoted to an interdisciplinary investigation, 
both in theory and experiment, of the hydromechanical, biophysical and some physio- 
logical aspects of muco-ciliary transport in the mammalian trachea, The tracheal muc: 
ciliary system is selected for its relative simplicity and its lack of other perplexing 
7 


factors such as the muscular contractile movement in oviduct. > biomedical 
cance of this research is relevant to certain diseases of and hazard to the respiratory 
function includ 


signili- 


tract that either are caused by or can result in abnormal muco-ciliary 
cystic fibrosis, chronic bronchitis, and adsorpti« f radioactive or | chemically- 


pollutants. 
The tracheal muco-ciliary system essentially consists of a mucous lay 
1] 


layer and epithelial cells. The epithelium is comprised of three main types of cells: 


cilia attached ciliated cells, mucus secreting goblet cells, and tl 


through a more refined classification further includes brush cells, serous cells and 


1e underlining basal cells, 


migratory cells (for a review see Jaffery and Reid [1]). Mucus, an aqueous suspension of 
the glycoproteins (2-4%) and salts (1-2% in water (95-97%) [2], is thought to have a tran- 
sition from a gel upper layer (5-10 um) to a sol lower layer (~6 pm) = 1e effects of 
ciliary activities increase toward the ciliated region [3-6]. The 

nd 


-elasti« shear-thinning medium under laboratory controlled condi 


muci has been obser 


to behave like a visco 
tions [7-9]. While a great deal of such experimental studies have been made, it 
stressed that adequate rheological properties of the mucus are generally n 
in vivo conditions with the strain and strain rate approaching the level 


states. On the other hand, the serous fluid is generally accepted t 
[8, 10] though its density and viscosity have not been well quantified. 
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The mu ‘ hown inside the trachea with and without a pressure gradi- 


ent. In the insert the « a density has been expanded to help illustrate the 
mucus, serous fluid, ar ilia interaction. 








Mucus transport : nd shear stress (b) for the rat trachea as a function of 
time, 


In a final experiment the was subjected to a pressure difference until the cilia 
action along the walls deformed > mucous plug sufficiently to cause the breaking of the 
plug and the formation of an annular coating on the inner walls of the trachea. When this 
ccurred SP immediately dropped to zero and the annular mucous blanket, which now 
more closely approximates the natural state of operation in the respiratory tract, moved 
at a velocity similar to the plug with AP QO as seen in Fig. 4(a) and (b), 
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TIME 


FIG. 3 Mucus transport rate (a) and shear stress (b) for the chicken trachea as a 
function of time. 





FIG. 4 Mucus (plug and annulus) transport rate (a) and shear stress (b) for the chicken 
trachea as a function of time. 
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where U_ is the maximum cilia tip velocity and K is a constant. Physically this signi- 
fies that fhe x-component of ity of the cilia away from the central station (x = 0) 


bution with 'wavelength' A= 2n7/K. We further 


17 


falls off according to this c 
denote 


phase angle 


ke a, Kx 


The boundary conditions + be prescribed 


a) , 


where 


IV ov 
+ =—) 


Ox 


oe 


rr _.? xr 


In addition, we have the no-slip « at the epithelial surface 


(r R 4 : {6 


16] vw) , (15) 
the last being the boundedness condition at the central axis, Here, conditions (8), (9) are 
in fact an approximation to the no-slip condition on the interfacial fluids in contact with the 
tips of the cilia patch for shallow depth of penetration. Condition (10) gives requirement 
on the continuity of velocity while nditions (11), (12) on the continuity of both normal and 
shear stresses across the part of interface which is free. These three conditions are the 
first approximation for small displacement of the free interface. 


In seeking the solution, it is « 


nvenient to adopt the Oberbeck potential represen- 
tation, by which u, v 


, and p can be written for the flow in the mucous plug 
(0 r R, -f 6< Ww) as 


I (r )] cos né , 
Oo n 


r )] sinng , 
n 


is the modified Besse! function of the first kind, 


(19) 


We see that P_ is proportional t 
o 

can be arbitrary, and U is equal t 

axis. 


the mean pressure gradient imposed on the plug, which 
the mean velocity of the mucous plug at the central 


For the flow in the serous layer, 


we have argued that the periodic component of zero 
mean can be neglected. 


Consequently we have for R r R +f, -w << OC, 
@) 
? 
+ A' log (r/R 0, (20) 
oO 
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l 
where (u/p')P. au" 


Substituting the expressions (16) - (22) in the boundary conditions (8) - (14), and 
subsequently using the approximation //R 1 and the asymptotic expansion of I (R_) for 
. . mn 


large Ro nkR = 2mnR/\ , we obtain the following dual Fourier series equations 


l 
> a cos né@ 
a n 


Thus, U and A can be determined from the dual equation (23), (24) and the solution to 
] 


the mucus flow is completed by the relations 


27) 


B -A_(4R_ - 3)/(4R_ + 3) 
n n n n 


The above dual Fourier-series equations can be immediately inverted in closed 
by Tranter's formula (Sneddon [31], p. 163). In particular, the constant coefficient 


found to be 


] 
=[P , (cos @)+ P (cos a@)] , 
c ve V 


where P (cos @) denotes the Legendre function of order v, v being given by (6). 

result provides the ratio of the mean velocity U of the plug at the central axis to the 
maximum cilia tip speed U.. The pressure gradient is seen to affect the velocity ratio I 
through an additive term ro as should be expected from the fact that this problem is 
linear. It therefore suffices only to discuss I , which is the velocity ratio at zero pres- 
sure gradient. 


The parameter @ 


As shown in Fig. 5 I depends on two parameters, @ and v. 
while v 


: . \e) 
is a measure of the size of the cilia patch in contact with the mucous layer, 
indicates the rate of decrease in longitudinal tip velocity of the off-center cilia from the 
If the patch size is small compared with the metachronal wavelength, 


maximum tip speed, 
range of moderate values of 1 From 


a@ is smalland | is not far less than 1 overa : 
oO 
our experimental measurement of the mean transport rate of an unobstructed mucus plug, 


we have estimated (see Eq. B) that 


U/T » 0.84 . 
c 
This suggests that if a@ n/3, which corresponds to the duration of effective stroke equal 
to 1/3 of the metachronal period, then v would lie in the range |! \ 2. While this 
estimate seems to be quite reasonable, critical comparison between theory and experiment 
will require some sophisticated optical systems to locate and track 


face, which has already proven to be a challenging feat. 


the muco-serous inter- 


From (13), (16) - (19) we obtain the first term in the axial shear stress evaluated 


R as 


l 
; 5 (31) 


xr ox 


This same expression when used in evaluating the mean 
- ? " 

5 - 21 N/m* for the chicken tracheal cilia. 

aP/ax 0) we would anticipate the 


which of course agrees with (A). 
cilia shear stress yields the estimate of 
Under physiologically normal operation (i.e. when 
actual shear stress to be somewhat less. 
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ABSTRACT 

A new model was developed to characterize the hydrodynamics of the interaction 

between spermatoza and cervical mucus. In this model the local environment 

about a length element of a spermatozoon in the mucus is represented as a 
consisting of a Newtonian fluid within a rigid, cylindrical no-slip 

boundary. Such a ceZZ models the mucous microstructure in the sense of 

“hydrodynamically equivalent” properties. A simple resistive force theory 

for sperm hydrodynamics was developed and applied to a human spermatozoon 

swimming within a mucous c: . The results included the prediction that 

the drag on the sperm head has a substantial effect on swimming speed. 

This prediction was tested experimentally, by analyzing high-speed cine 

films of normal and headless human spermatozoa swimming in mucus. The 

theoretical result was confirmed, and a size of the order of 2 um 

was suggested. 


INTRODUCTION 


Immediately following insemination, the spermatozoa of many mammals must negotiate a 
mucus-filled cervical canal as they begin their migration toward the egg (1,2). The cervix 
thus performs a crucial role in mammalian sperm transport, acting during the periovulatory 
period to admit and possibly to store normal viable spermatozoa, while simultaneously 
blocking the entrance of morphologically abnormal spermatozoa and pathogenic material. 
During non-ovulatory phases of the female reproductive cycle, and during pregnancy, no 
spermatozoa or other material are admitted. Much of this regulatory capacity of the cervix 
can be attributed to its mucous secretion. The cervical mucus, like the secretions of 
other epithelial tissues, is a hydrogel. Ninety to 99.5% of the mucus consists of the 
aqueous phase, which is similar in composition to blood serum and contains soluble proteins 
and inorganic salts. The solid phase is comprised of a carbohydrate-rich glycoprotein 
(the mucin) of high molecular weight. The three dimensional organization and integrity of 
the mucous microstructure depend upon interacting electro-chemico-physical properties of 
the mucin and mucous plasma (3). 


The migration in vivo of spermatozoa into and through cervical mucus is a transport 
process whose mechanics depend upon: the active undulations of the sperm flagellum; the 
morphological and possibly electrical properties of the sperm body; the micro and macro- 
physical characteristics of the cervical mucus as time and space-dependent entities; the 
geometry of the cervical canal, the pattern of mucosal folding and distribution of cilia; 
and the activity of the female reproductive tract, including the ciliary beat (where 
present) and muscular contractions of the viscera. Hydrodynamic interactions among these 
factors, with physical and chemical modulation at the molecular level, determine the 
swimming trajectories of spermatozoa during passage through the cervix. It is presently 
believed that active swimming by spermatozoa is a necessary if not sufficient component 
of the mechanics of sperm transport through the cervix (2). Attempts at describing this 
element of the sperm-mucus interaction should characterize both the sperm and the mucus 
in terms of parameters that relate to first principles of biomechanics and that are 
amenable to experimental measurement. 


There exist a small number of theoretical studies directly address 
mental hydrodynamics of the sperm-mucus interaction (4,5,6). These pre 


d to the funda- 
iminary treatments 
modeled the sperm body as a two-dimensional, small amplitude wavy sheet of infinite length. 
The representation of mucous properties included the presence of adjacent undulating sheet 
as a model of the mucin macromolecules. Not surprisingly, these studies indicated that the 
nearby presence of the mucous microstructure could significantly influence the flagellar 
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The pooled data are illustrated in Figure 2, in which the dimensionless swimming speed or 
"kinetic efficiency” V,/fL has been plotted vs. beat amplitude. Here L is the length of the 
sperm body, which is 50 um for a human spermatozoon. The straight lines in Figure 2 are 
linear regression lines fitted through the origin. The slopes of these lines are 1.17 for 
the normal sperm, and 2.49 for the headless sperm, and they differ significantly (P < 0.01). 
The legitimacy of this comparison procedure was established by determining that at a 5% 
significance level the unconstrained regression lines were linear and had intercepts not 
Significantly different from zero These data indicate that a headless spermatozoon tends 
to swim approximately twice as fast in the mucus as a normal sperm exhibiting identical 

beat frequency and amplitude (viz. beat shape). 


RETICAL DEVELOPMENTS 


As an initial new approach for describing the hydrodynamics of the sperm-mucus inter- 
action, a simple tt force t) was developed (14). The mucous microstructure 
experienced locally by length elemer among a sperm body was represented by a rigid cylin- 
drical no-slip boundary. This model of the hydrodynamic influence of the microstructure was 
motivated by our contemporary, qualitative understanding of its geometrical properties. The 
mucin network in the mucous native gel state is comprised of filamentous macromolecules 
whose length and diameter are smal] compared to a sperm flagellum. These macromolecules 
are organized such that the interstices among them vary in shape and orientation. A 
typical interstitial size is also all compared to the length of a sperm body. Thus, a 
spermatozoon swims in very close proximity and interacts hydrodynamically with numbers of 
the mucin macromolecules. The local force on each length element of the sperm body is 
assumed to depend upon the velocity of that element only. Such an approximation is rela- 
tively accurate for a structureless Newtonian fluid of infinite extent (15,16,17). The 

icin macromolecules may, in fact, shield local elements of the sperm body from direct 
hydrodynamic interaction, therefore increasing the rationale for the resistive force 
approach. 


The following resi: yrce efficients were employed, cf. (18): 


Here C; and Cy are respectively longitudinal and normal force coefficients for the 
flagelTum, and Cy is the force coefficient for the sperm head, assumed spherical; n is the 
viscosity of the surrounding fluid, assumed Newtonian; ro and ay are respectively the radii 
of the flagellum and head; and h i he separation distance between the centerline of the 
sperm body and the outer boundary se of these expressions is predicated on the assumption 
that the sperm swims undirectionally along the centerline of the cylindrical cavity. It 
follows that a simple analytic expression can be derived for the propulsive velocity or 

4) 


kinetic efficiency, 4 viz. 


. 


is an effective j 2S measured along the flagellum; and 8 is the mean-square 
value of the cosine of loc e angle of the flagellum with respect to the axis of 
propulsion. 


Tr 


In order to model f beat of human speramtozoa in cervical mucus, the 


waveform depicted in Fi i. ised. The sperm head, if present, was represented 

by a sphere of radius 1.25 um. é agellum was taken as a circular cylinder of length 
45 ym and radius 0.25 um. The proximal 10 wm segment of the flagellum, which undergoes 
little bending, was modeled as a rigid circular cylinder. The distal 35 um of the flagellum 
was then taken to propagate traveling sine waves. Three combinations of amplitude (b) and 
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wavelength (A) were applied for these waves: b = 2 um, A = 25 um; b = 4 um, A = 25 um; 


and b = 5 wm, A = 25 um. Here A is the wavelength as measured along the axis of propaga- 
tion. 


Figure 4 displays the results of the application of these waveforms to Eq. 4 for 
spermatozoa with or without heads. The theoretical model predicts that for cell sizes of 
the order of 2 um, headless spermatozoa would swim approximately twice as fast as normal 
spermatozoa with heads and exhibiting the same beat. 





FIG. 3 
The model beat shape used in the theoretical calculations of a human sperm swimming 
in cervical mucus. 


FIG. 4 
Results of the theoretical model of a human spermatozoon swimming in cervical 
mucus. The symbols are: V, = propulsive swimming speed; f flagellar beat 
frequency; b = beat amp1ituge; k = 2n/(beat wavelength); L = length of sperm 


body (50 um); and h =cell size. The broken lines are for headless spermatozoa. 
The solid lines are for normal sperm with heads. 


DISCUSSION 


The swimming of a spermatozoon within cervical mucus is a phenomenon in which hydro- 
dynamic interactions between the sperm and the mucous microstructure are strong. This 
can be appreciated by reference to the theoretical model developed here with, for example, 
a oc size h = 2 um, a beat amplitude b 4 um, and a wavelength ) 25 um. In such a 
case, the model predicts that a normal sperm (with a head) would have a kinetic efficiency 
(dimensionless swimming speed) approximately 50% higher in the mucus than while swimming 
with the same beat in a simple Newtonian fluid of infinite extent. This distinction arises 
from the amplification within the mucus of the local forces acting on the sperm body. For 
example, the normal resistive force coefficient Cy increases threefold in the mucous 
as compared with unbounded Newtonian fluid of the same local viscosity. The consequent 
increase in local bending moments acting along the flagellum undoubtedly contributes to the 
altered beat shape of mammalian sperm in the mucus as compared to simple Newtonian media 
of low viscosity, see e.g. (19,20,21). The increased drag on the sperm head, resulting in 
an altered “boundary condition" for the flagellar contraction mechanism is also likely to 
be a contributing factor (22). For example, the drag coefficient on the sperm head 
increases fourfold in the mucous ce/2 as compared to an unbounded Newtonian fluid of the 
same local viscosity. Thus, the tradeoff between the increased drag on the sperm head and 
the increased forces acting upon and generated by the flagellum may be directly involved 
in modulating both the beat shapes and swimming speeds of spermatozoa in cervical mucus. 
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LETTER TO THE EDITOR-IN-CHIEF 





COMMENTS ON "A TWO-FLUID MODEL FOR BLOOD FLOW 
THROUGH SMALL DIAMETER TUBES" 


YVES NUBAR 
14 Horatio Street, New York, N.Y. 10014 


May I refer to a recent paper by Chaturani and Upadhya (1), a theore- 
tical study of the flow of blood viewed as a composite fluid with a couple 
stress fluid core and a peripheral shell of plasma. This is a valuable work, 
well organized and discussed with care. One is however somewhat jarred by 

the unconcern and casualness with which the authors set the flow velocity 

at the boundaries equal to zero. This unquestioning espousal and use of the 
handy "no-slip" hypothesis, the simplest possible of all boundary conditions, 
is not uncommon among students of fluid flow, due perhaps to ingrained habit, 
prevailing attitudes, or the path of least resistance. The fact is, however, 
that a number of studies of blood flow, both theoretical and experimental, 
have suggested, some by implication, others outright, the likely presence of 
slip (or a velocity discontinuity) at the flow boundaries (or in their imme- 
diate vicinity) - thus making the "no-slip" hypothesis one of uncertain 
validity. 


Among these studies - some of them discussing slip in some detail and 
at some length - may be cited those by L.E. Bayliss (2,3), L. Bennett (4), 
E.H. Bloch (5), G. Bugliarello and J.W. Hayden (6), A.L. Copley (7), I. Isen- 
berg (8), A.L. Jones (9), A.B. Metzner (10), Y. Nubar (11,12), G.W. Scott 
Blair (13), V. Vand (14). H. Bateman (15) and S. Goldstein (16) give summaries 
of the history of slip in addition to discussing boundary conditions in general. 


If slip at the wall is indeed present in the flow of blood, as appears 
probable, ignoring it as a boundary condition in theoretical flow studies 
would be inaccurate. This applies equally in the case of experimental, visco- 
metric studies of blood flow in which slip in the viscometers would affect | 
the readings materially, unless it is detected and corrected for - something 
rarely, if ever, attempted. This is doubly inopportune, since these readings 
are normally used as criteria in testing and evaluating theoretical studies 
of blood flow. 
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ALAN CHADBURN BURTON (1904 - 1979) 


Alan Burton, M.B.E., B.Sc., M.A., Ph.D., LL.D., D.Sc., 
F.LR.S.C., F.R.S.A.(U.K.), Professor-Emeritus at the University of 
Western Ontario, and one of Canada's distinguished biophysicists, 
collapsed and died June 27, 1979 on his way home from work. For 
those of us who knew Alan Burton over the years it is a time to 
reflect on him as a scientist and teacher and on what his special 
contributions were. 


He was born in London, England, grew up in a large middle 
class family, and through two wars and the Great Depression, was 
exposed to a spectrum of social changes and economic circumstances. 
Variety was important to him in his life and scientific work. He 
received First Class Honors in Physics and Mathematics from the 
University of London, and came to Canada to do postgraduate work 
in Physics at the University of Toronto. He received his M.A. for 
research on ‘The absorption spectra of the major planets', and his 
Ph.D. for his research ‘Heating of electrolytes by high frequency 
currents'. His entry into the world of biology and the start of 
his career as a biophysicist, was, in his words ‘more or less 





Portrai of Alan Burton, on the occasion of his retire- 
ment, drawn by Dr. Alfred Jay. Dr. Jay wi graduate 
student of Dr. Burton. 





accidental’. He began by explaining to physiologists and physicians 
the heating of body tissues by shortwave radiation. Alan Burton was 
quickly attracted and committed to the physical challenges of 
biology; he trained with Professor Martin at the University of 
Rochester to work on heat exchange in man and animals, and later 
with Professor Bazett at the Medical School of the University of 
Pennsylvania to study peripheral blood flow, followed by a research 
fellowship for four years at the Johnson Foundation for Medical 
Physics at the University of Pennsylvania, working with Dr. Detlev 
Bronk on human temperature regulation. At the outbreak of World 

War Il he returned to Canada as a Research Associate of the National 
Research Council to work at the Banting Institute, Toronto. For his 
war researches in aviation medicine and on protective clothing he 
was awarded the M.B.E. 


After the War, Alan Burton was invited to join the Depart- 
ment of Medical Research at the University of Western Ontario. He 
came, and three years later, was appointed Professor and Head of the 
new Department of Biophysi in the medical school. He was strongly 
committed to postgraduate education; the University, with his 
leadership, established the first Canadian postgraduate program in 
biophysics. Over the years he personally trained 14 M.Sc. and 22 
Ph.D. students. He worked to establish Canada's first Honors 
Biophysics program at the iversity which began in 1966. 


Alan Burton earned an international reputation as a 
scientist and received many awards in recognition of his scientific 
contributions. He was elected a Fellow of the Royal Society of 
Canada (1950), President of the American Physiological Society 
(1956), President of the Canadian Physiological Society (1959) and 
President of the Biophysical Society (1966). He served as Chairman 
of the Federation of the American Societies for Experimental Biology 
(1957-58) and of the Canadian Federation of Biological Societies 
(1963). He received sairdner Foundation Award (1961) and the 
degree (Honoris Causa) of LL.D. from the University of Alberta (1964) 
and D.Sc. from the University of Western Ontario (1974). The major 
areas of his published work were temperature regulation and blood 
flow, the equilibrium mechanics of the blood vessel wall, the mecha- 
nics and geometry of erythrocytes, the deformability of membranes and 
permeability, and in recent years, intercellular communication and 
the theories and epidemiology of cancer. He was perceptive about 
basic mechanisms in science and enjoyed writing. His book, published 
with Edholm in 1955, ''Man in a Cold Environment" is regarded as a 
classic. His second book ''The Physiology and Biophysics of the 
Circulation'', published in 1965 (and revised in 1972) has been widely 
used by medical students and was translated into French, German, 
Italian, Japanese and Spanish. His third book, which draws together 
the thrust of his research on cancer, will be published late in 1980, 
“Understanding Human Cancer: The Physiological and Biophysical 
Point of View''. 


During the years since his retirement in 1970 he continued 
to teach. He vigorously pursued his new field of research on cancer 





and argued his theories with colleagues and lectured about his ideas 
at other universities. He was uncomplaining about his declining 
health, and his deteriorating eyesight. 


Association with Alan Burton over the years leaves an 
accumulation of memories of a rare man, a man whose scientific 
perception commanded international respect, and yet whose humility 
made him accessible to every level of student. He would give 
willingly of his time for a frustrated student and never judged 
anyone's efforts unworthy of his attention. Because he was always 
an avid learner, he was a fine teacher who welcomed questions as an 
opportunity for a fresh examination of the material and felt it his 
responsibility to bring the students to a comprehensive understanding 
of his subject. In his department, coffee breaks were obligatory 
periods of rigorous training for all aspiring scientists. He would 
range about the group requesting up-dates on experimental progress, 
assessments of lectures attended, input on problems encountered by 
any individual, and frequently would present aspects of papers he 
was reviewing as targets for departmental brainstorming. On these 
occasions the professor would join right in and exercise his own 
formidable and wide-ranging talents -- sometimes forcing students 
and faculty alike to return to fundamental concepts to find the 
significance of amassed data, making an original attack on old 
material with a deceptively simple question, seeing relationships 
in apparently disassociated studies, and allowing his humorous 
side to create a party atmosphere so all could respond to his 
challenges in the spirit of shared fun. He would often draw 
hilarious comparisons, invent rhyming couplets, plunge into an 
impromptu demonstration with any available materials such as 
cigarettes, matches, or coffee cups (with messy if not disastrous 
results), and thus infect everyone with his joy in science. 


The study of science so enriched Alan Burton's life 
that he was tireless in initiating others in this fraternity. 
From science he derived challenge, purpose, values, adventure, 
worth, and profound joy. He lamented the modern trend to special- 
ization as dessicating the pursuit of science. He welcomed 
variety and change and saw a growing need to reach across the 
isolation of specialties and synthesize facts into true knowledge. 
His most treasured quotation was Francis Bacon's description of a 
scientific mind:? 


In 1975, Alan Burton said, "It has been my aspiration to live up 





to Bacon's description, which explains, perhaps, how anyone could 

so spread his efforts on such a variety of topics in science.''” 

He considered himself a most fortunate man to have been able to 

live his life as he did. And we have been fortunate to share in it. 


Peter B. Canham, Ph.D., 
Professor of Biophysics, 
University of Western Ontario. 


‘Selected Writings on the Interpretation of Nature. Francis Bacon. 
Random House, New York. 1955. p. I51. 


“Variety -- the spice of science as well as of life. The disadvant- 
ages of specialization. A.C. Burton. Ann. Rev. Physiol. 37, 1975. 
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FLOW OF BLOOD IN TRANSFUSION SYSTEMS 


Palosaari* and K. Linko** 


* Department of Chemical Technology, Lappeenranta University 
of Technology, Box 20, Lappeenranta 85, Finland 


** Department of Anaesthesia, Helsinki University Central 
Hospital, Hoartmaninkatu 4, 00290 Helsinki 29, Finland 


ABSTRACT 
The flow resistances caused by commercial elements used in blood transfusion 
systems were measured. These elements are comprised of microfilters, blood warmers, 
transfusion sets, and a venous cannula. In the absence of any occlusion in the 
filter it was found that the flow resistances coused by microfilters were low and 
without practical importance. Two warming devices were tested. The Fenwall® “dry- 
heat" warmer was found to create a high flow resistance whereas the Portex® Coil 
created only one third of the resistance of the Fenwal "“dry-heat" warmer. Calcula- 
tions from the experiments show that it is possible to substantially reduce the 
overall flow resistance created by the transfusion set and warmer. This may have 
practical importance although in practice the limiting factors to the flow may still 
be occlusion in the microfilter as well as flow resistance in the veins of the 
patient. In this work ao simple method wos developed enabling the flow resistanc 
measurements to be made by the use of water or any Newtonian liquid. These values 
can be used to estimate the flow resistances to blood. 


INTRODUCT ION 


Previous results on the flow rates of aggregate-free blood [| 1, 2 }] suggest that the flow 
resistance caused by in-line blood warmers limits their usefulness. Since these devices ar 
already in general use it is important to study their flow properties in more detoil. It is 
also necessary to acquire more information concerning the resistances of other elements of 
the transfusion system and to discuss blood flow in order to find out ways to minimize th 
overall flow resistance. 


According to experimental evidence by Walburn and Schneck [ | yeneral power law, i.« 
? 


’ 
the Ostwald-de Waele equation is valid for blood for shear rates sve 23 reciprocal seconds 
9 


n r . 
where K is the flow consistency index, Pas ;1 is the shear stress, Pa n is the flow 


behoviour index, dimensionless ; and dv/dy is the shear rote, i.e. é ite of chang 
velocity of the liquid with the distance from the solid wall. 


Often the concept of apparent viscosity , r , is being used. is fined as follows 


app 


Walburn and Schneck [ 3 |] developed various experimental expressions for the quantities K and 
n of equation (1) for blood. For the accuracy sufficient in the present work, the following 
expressions can be used, for blood at a temperature of 37°C. 
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the tube diamete e rate of flow, and po is the density of the fluid. 
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ann , : - 
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Because blood is not a Newtonian liquid, Fig. 4, as such, cannot be used to estimate the 
blood flow resistances of the various elements. According to equation (5), the flow 
behaviour index, n, appears in the exponent of the diameter of the flow conduit, which is 
not known. To overcome this difficulty it is necessary to determine the flow resistance of 
the element for two liquids of different viscosities. Subsequently, the equivalent lengtl 


and equivalent diameter of the flow channel can be calculated from equation (5). This 


method was adopted and the length and the internal diameter of a venou: mnula was 
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calculated from the pressure drop measurements with glycerol-water and blood. Unfortunately, 
this method proved to be rather sensitive to measurement errors, and erroneous values: wer 
obtained for the ength and diameter, i.e. slight inaccuracies in measuren reate gr 


inaccuracies in length and diameter. 
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BIORHEOLOGICAL ASPECTS OF BLOOD FLOW THROUGH ARTERY WITH 
MILD STENOSIS : EFFECTS OF PERIPHERAL LAYER 


J.B. Shukla, S.P. Gupta and 2.S. Parihar 


Department of Mathematics ,Indim Institute of Technology 
Kanpur-208016, INDIA. 


The blood flow through an artery with mild stenosis has been analysed 

by characterising it as a power-law fluid. It has been shown that the 
resistance to flow and the wall shear stress, as comared to non-stenosis 
case, decrease as the consistency of peripheral layer fluid decreases 

and this decrease is ehanced due to the psewioplastic nature of the 
blood, 


INTRODUCTION 


Several workers have stulied the flow characteristics of blood in an artery with mild 
stenosis by considering it as a Newtonian fluid [ Young (1,2), Morrester and Young (3)]|. 
However, it has been observed that blood, being a suspension of cells, behaves as a non 
Newtonian fluid at low shear rates in tubes of smaller diameters | Hershey et al. (4), 
Huckaba et al. (5), Charm and Kurlani (6), Whitmore (7), Han and Barnett (8)|. it has been 
reported that the suspesion of red cells in serum, the behaviour of which is similar to the 
whole blood for shear rates more tha 1 sec~!, obey power-law mdel [ Whitmore (7), Han and 
Barnett (8)] - Han and Barnett (8) pointed out that the power law exponent for blood serum 
may vary between .2 to .9 depending upon the pathological conditions, while Whitmre (7) 
quoted this value between .68 - .80. The exponent for plasma has bem reported to be .95 - 
-995. The consistency of the suspension of red cells in serum is found to be ten times the 
viscosity of the suspending medium [ Whitmore (7), Pe 73 | ° 


Keeping in view the non-Newtonian power-law behaviour of blood, its flow characteristics 
in an artery with mild stenosis have been investigated by Shukla et al. (9). However, they 
did not take into accoumt the effects of rajiial distribution of cells and the existence of the 
peripheral plasma layer near the wall whose thickness may be of the order 40m in tubes of 
diameter 1 to 6 mm | Hershey et al. (10),Lih (11), Middleman (12)] - In this p@er, an attempt 
has been maie to study these effects on the characteristics of the blood flow in an artery 
with mild stenosis. 


ANALYSIS 


Consider an axially symmetric flow of blood (fluid) through an artery (circular tube)with 

mild stenosis (constriction) as shown in Fig. 1. The blood flow is assumed to be characteri_ 
sed by the power-law model whose consistency varies symmetrically along the radial direction. 
As pointed out by Caro et al. (13, p. 111), the time average value of the shear would have 
Similar spatial distribution as associated with stealy flow, hence, the flow studied here is 
assumed to be laminar and steady. 

R( 2) Os 2% “o 

5 = 1 = Se [t + cos F(a - a 22)]; a fz Clr 

° ° ° 


=1 ; elsewhere 


(1) 


where R(z) is the rafius of the artery in the stenotic region, R is the constant radius else- 


where, L, is the length of the stenosis and 6, is its maximum height (4, << Ris Pig. 1). 
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7eometry of arterial stenosis 


equation governing the 


is | Shukla et al. 9) } 
38 = 

am (97 {) « 

re) 


power-law fluid is given 


(3) 
w is the axial velocity, is the fluid pressure, mr) is the consistency function ani no is 
the power-law exponent, assumed to be a constant and may be taken as the mean of the flow 
behaviour indices over the radius, 
olving equation md using the boundary conditions 


the axial velocity w is 


The volumetric flow rate, Y, is 


which on 


dr. 
/ >! n 
° Sn(r 
L J 
It can be seen from equation of continuity that the flow flu 
section of the tube. 


2 is constant across any 
Rearranging equation { 





Integrating equation 
bave 


dz 
. n* 
{1(z)} 
The resistance to flow, A, is defined as 


p - D 


oe *o 


which o0 using equation (10 


In the case of no stenosis 6, O or Riz) = R ) and with constant consistency of 


4 


fluid (m= a)» say), the resistance to flow, wy? is calculated from equation (12) as 


»y 


1 





dz. 
30 
{1(z)} 


The shearing stress, I ., on the wall of the tube is given by 


Own | 
=) : 
+r=R( z) 


which on using equations (5) « 8 


gives the shearing stress, t» at the maximum height of 
the stenosis (i.e. at z /2), a8 


The wall shearing stress, T., 


calculated from equation (15) as 


in the case of no stenosis with the constant consistency m, is 


‘< 
°F (17) 


Prom equations (16) and (17), we have 


z +1 
pr a 
n Bp °s, 0 
/ 


= ( — i 7 
5 Ty Gne1? Weg m, mer |. 


The resistance to flow and the wall shear stress at the maximum height of the stenosis 
can be stuwiied from equations (14) and (18) respectively for any general consistency function 
mr). 





EFFECTS OF PERIPHERAL LAYER 





as pointed out earlier, the existence of cell depleted layer near the arterial wall has 
been suggested by many workers| Hershey et al. (10), Middleman (12)] - To see the effects 





of consistency of such a peripheral layer fluid on the flow behaviour, the function m(r) may 


be assumed as follows 


(19) 


; 8 r< R 


< 
2 l=- _ 
2 


where m,, ™, are the consistencies of the fluids of the central and the peripheral layers 
* 


respectively. It may be noted that the value of the power-law exponent, n, in this case is 
taken to be the mean of the flow behaviour indices of the central and peripheral layers. 


r 


Peripheral layer 








FIG. 2 


;eometry of arterial stenosis with peripheral layer 


The function R,\2 represaits the geometry of the interface ami may be given by (see 
Fig. 2 
R,\ 2) ‘ _- L 
1 11+ cos = (2-4 ~ 2) ; a<cs CL + 
R L L 2 = 
° ° 

(20) 
= ; @leewhere 


where q is the ratio of the central layer raiius to the tube radius in the unobstructed 
region and 6. is the maximm bulging of the interface due to stq@osis which is to be 


.e) 


determined as follows 


Using the consistency fumetion (19), the fluid velocities w,, ¥ of the central and the 
peripheral layers respectively have been calculated from equation (5) to obtain the 
corresponding fluxes % and > as 


. 


f ear, dr 


The total flux, Q is 
Q = So + % 





+ 3 a i "7 
oa €4 og | (23) 


This flux Q cam also be obtained from equation (7) on using the consistency function (19). 


It can be se@m by integrating the equation of continuity appropriately that not only Q 
but & and Q are also constants separately. 


Now, integrating equations (21), (22) ami (23) and remembering that the pressure drop, 
(Py - P, across the tube length is the same in each case, we get 


—- 


gi/n 3n+1 ( 











(29) 


Since Q= 4 + +9 from equations (24), 
to determine Og» 


5n+1 nti 





It can be seen by direct substitution that 


Ry = aR 
Satisfies equation (30), which on using equations (1) ami (20) gives 
O4 “Gdg° 
Using the consistency function (19) in equation (9) and making use of the relation (32) 


in equations (14), (18), the expressions for 4, A Can be written in a simplified form as 
follows 
A= { 


gi/n 





bettas” 











It may be noted that for general value of n, the integral (35) cannot be evaluated 
analytically. But whem 3n+1 is an integer, this integral cm be evaluated by using the method 
of calculus of residues. mn such a case, A, from equation (33), can be written as 





It is seen that equations (34 and (56) reduce to the case of Shukla et 
1.0 and to the case of Young (1) fora 1.0, 8 = 1.0. 


RESULTS AND DISCUSSION 





To see the effects of various parameters on /* Tt,» the following values have been used 


Whitmore (7),Han and Barnett (8), Hershey et al. (10), Middleman (12)] $ 


a ).667, 


» = Oel, 0.3, 


“o 


in Fi 3 ami 4. From these figures, it is noted that 


Using these data the values of x from equation (33) or (36) have been calculated 
Oo 


numerically and plotted against 6 R 


i increases as the stenosis size 6. s increases for fixed § and n, It is also 


° 
noted from Fig. 3 that \ decreases as 8 or decreases for fixed stenosis size. 


c 


G. 3 


Variation of } with 6 /R, for 
different n and 8 














Variation of , A with 6, R, for different n and & 


The values of .. from equation (34) have been plotted in Fig. 5 for various 


B and Oo. R,° It is seen from this figure that 7, also increases as the stenosis 
increases but decreases as 6 orn decreases. 


The flow characteristics of blood flow through an artery with mild stenosis have been 
investigated by modelling the blood as a non-Newtonian power-law fluid with variable consist- 
ency, By assuming the existence of a peripheral layer with a fluid of smaller consistency, it 
has been shown that the resistence to flow and the wall shear stress, as compared to non- 
stenosis case, (though increase with stenosis size) decrease as the consistency of the 
peripheral layer fluid decreases md this decrease is enhanced due to the pseuio-pla 
nature of the fluid. 





Thus, it appears that nature has provided a self controlling mechanism by associating 
favourable conditions to the system so that it may function smoothly as far as possible even 
under diseased conditions, 
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The time dependent theory of Arminsk Weinbaum 

the vesiculi transport of acromolecules across 

modified to include the effects of steric hindrance of ve at- 
tached at the plasmalemma membra f the cells. The spatial varia- 
tion of the diffusion coefficient near the plasmalemma membranes of 
an endothelial cell layer is determined experimentally by measuring 
velocities of spheres settling under gravity perpendicular to 

gle plate fitted with an array of attached spheres th mode] 
trastructural features of plasmalemma membrane and 

vesicles. The spatial v ation of the diffusion coefficient in the 
central part of the cell is determined by the recent hydrodynami 


theory of Ganatos, Weinbaum and Pfeffer (2) for the creeping motion 


of a sphere between two plane parallel walls. This improved descrip- 
tion of the vesicle dynamics is used to generate vesicle concentra- 
tion profiles using the theoretical formulation of Arminski et al. 
(1). A comparison of the present solutions with the time dependent 
tracer data of Casley-Smith and Chin (3) shows substantially improved 
agreement over previous models in which the effect of the steric hin- 
drance of attached vesicles is neglected. 


i. INTRODUCTION 


Electron microscopic studies (4,5) and perfusion st 
shown that the endothelial cell layer is the rate limiting 
transport of macromolecules across the artery wall. 
scription by Palade (9), there has been increasing evidence 
transendothelial transport of macromolecules, especially those 
40A in diameter, cross the cell by vesicular transport. 

With some deviation for vascular heterogeneity, electron dense markers 
larger th 40k are seen to traverse endothelial cells via f sma la vesi- 
cles rather than through the intercellular clefts (5). Ves 
cal, 700%, membrane bound intracellular bodies. When unimpeded by cytoplas- 
mic structure, unattached (free) vesicles unde: a diffusional migratior 
across the cell which is strongly influenced by a hydrodynamic d molecular 
force interaction with the plasmalemma membranes of the endothel 
The time dependent tracer studies of Casley-Smith and Chin c 
Simionescu, Simionescu and Palade (5) indicate that vesicles are 
long enough at both the luminal and tissue fronts for the marker 
horseradish peroxidase and exogenous myoglobin to come to concen 
librium with the surrounding luminal 
Current evidence indicates that larger 
do encounter molecular sieving 
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Schematic illustration of experimental apparatus. ( 
tube support removed, ruler, mirror and movie camera. 
tank. 


and attached spheres were designed to model the ultrastruc 
plasmalemma membrane and its attached vesicles. The spher 
in. corresponded to a vesicle diameter of 7OOA. All other dimensions were 
scaled according to these figures, including the length of he connecting rod 
between the plate and attached spheres so that a vesicle neck was simulated. 

Three different plates were used each with a different center-to-center 
spacing between attached spheres. The at 1 center-to-center 
distance, y, and vesicle number density per um, ©, isshown in Fig. 4. A 
hexagonal pattern was chosen for the array after carefully viewing the 
freeze-cleaved preparations of Simionescu, Simionescu and Palade (18) showin: 
broken vesicle necks (see for example, Fig. 2 in that paper). 

The tank was placed on a fatigue mat in order to damp out vibrations. 
Care was taken to insure that the plate was parallel to the surface of the 
fluid. The experimental set-up was maintained at constant temperature so 
that the viscosity and density of the fluid were nearly uniform. Temperature 
and viscosity measurements were made at the beginning and end of each series 
of runs to insure that there were no changes in viscosity during the runs and 
no temperature gradients in the tank. 
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touched an attached sphere; thus only hydrodynamic effects were considered. 
4. RESULTS AND DISCUSSION 
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CHANGES IN SURFACE AREA AND VOLUME 
MEASURED BY MICROPIPETTE ASPIRATION 
FOR ERYTHROCYTES AGEING IN VIVO, 


G.B. NASH and S.J. WYARD 
Department of Physics, Guy's Hospital Medical School 
London Bridge SE1 9RT England. 


INTRODUCTION 


During in vivo ageing, erythrocytes decrease in volume and lipid content whilst their mean cell 
haemoglobin concentration (MCHC) and density increase (1, 2, 3). At the same time their enzymatic 
and metabolic activities lessen (4), The variation in density has made it possible for populations 
of red cells of differing age to be separated centrifugally, either using density gradients (2, 5) 
or by sedimenting whole blood (3, 6) and the success of the separation has been tested by cohort 
labelling of cells with radioactive iron and glycine. 


The dominant process by which senescent red cells are removed from the circulation is not 
certainly known, It has been suggested that a decrease in the deformability of aged cells may 
impair their ability to pass through the microcirculation, especially of the spleen, hence 
facilitating their removal by that organ (7, 8). In general, erythrocyte deformability depends on 
the elastic properties of the plasma membrane, on the fluidity of the cell contents and on the cell 
surface area/volume ratio (and hence shape) (9). It has been shown that aged cells require greater 
pressure to force them through narrow pipettes ( < 3 um internal diameter) (10), and that the 
rate at which they initially enter similar micropipettes is reduced (11), It has also been 
demonstrated that membrane fluidity, assessed by ESR methods, is lower in older cells (12). 

In addition it has recently been demonstrated taat the increase in MCHC with age is linked to a 
disproportionately large rise in the internal viscosity of the cells (13) which would decrease the 
rate at which they can be deformed, 


With respect to changes in size and shape, the decrease in cell volume which accompanies 
ageing has been well documented (2, 14). It has been suggested that erythrocyte fragmentation 
(10, 15) occurs which would be accompanied by a loss in membrane area and surface area/volume 
ratio. This would result in a decrease in the minimum cylindrical diameter which an erythrocyte 
could attain and might lead to mechanical lysis if a blood vessel of smaller diameter was 
encountered (7), In any case, sphering of red cells has been shown to decrease their ease of 
deformation (9) and it has been reported that separated aged cells appear to lose their biconcavity 


(12), 


In this study we have carried out direct measurements of individual red cell surface 
volume, in order to see if quantitative shape changes occur during ageing. A micropipette 
aspiration technique similar to that of Jay & Canham (16) was used to measure these parameters 
for erythrocytes confined to a cylindrically symmetrical geometry. Both human and rabbit blood 


were fractionated according to density to provide populations of cells of differing age. 


KEYWORDS: Erythrocyte Ageing, Erythrocyte Volume, Erythrocytes 
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Measurements on Separated Cells 


Aliquots from the cell suspension re placed in a glass chamber on a microscope stage. The 
chamber was closed from above with a erslip, and a micropipette, mounted on a motorised micro- 
manipulator (Micro Instruments Ltd., xford, England), was introduced via a channel in its side. 
Cells and pipette were observed using a X40, long working distance, water immersion lens (Carl 
Zeiss Ltd., London). Photomicrograph re recorded on 35mm film, The pipette was connected to 
a liquid filled, spring loaded reservoir which was moved vertically by a micrometer barrel. The 
pressure head was read directly from the micrometer scale after zero pressure had been adjusted by 
observing the behaviour of cells adjacent to the pipette tip. Micropipette internal diameter, approx- 


imately 1.8-2,.0um, was measured by h yn photography using a X40 objective lens before the 


pipette was filled with Isoton (Coulter Electronics Ltd., Harpenden, England). 


The correct pressure for deforming t erythrocytes was obtained by gradually increasing the 
suction for a trapped cell. Initially cell nd to fold as they enter the pipette but at a negative 
pressure of about lcmH,O (depending on pipette diameter) a smooth outline with a spherical outer 





and cylindrical inner portion is obtained. Having attained the correct suction, the reservoir was set 
at this level and cells were photographed a few seconds after entrapment. Subsequently they were 
ejected by raising the reservoir against its spring. The effects of further increasing the suction 
pressure and of prolonged cell attachment were tested to see whether these would alter the cells’ 
behaviour. Later the developed film negatives were projected and the diameters of the spherical 
portions and the lengths of the cylindrical portions of the cells were measured, From the values 
obtained the surface area, volume and surface area/volume ratio were calculated for each cell 
(assuming the tips of the cylindrical portions to be hemispherical). In addition undeformed cel] 
diameter was measured for erythrocytes settled on the bottom of the chamber. 


A Coulter Counter Model ZF with Channelyser and volume distribution plotter (Coulter Elec- 
tronics Ltd.) was also used to measure cell volume and to count the cells in the fractionated samples. 
These were diluted with KRBG and compared to freshly suspended whole blood, Volume calibration 
was carried out using latex spheres (median diameter 4. 9um) and thus to calculate absolute cell 
volume a shape factor was required for the cells (18, see also Results). As the sphere diameter 
quoted by the manufacturers was the median value, we decided also to use the median value for the 
red cell volume. Finally reticulocyte counts were made microscopically after staining concentrated 


cell suspensions with brilliant cressyl blue, by counting approximately one thousand cells. 
RESULTS 


Blood samples from each of three rabbits and three humans were fractionated according to 
density on two occasions, The more dense cells had markedly smaller volumes and reticulocytes 
were concentrated in the uppermost layers (see Table 1), The separation procedure did not alter 
overall median cell volume compared to freshly suspended cells (see Table 2) and the fractionated 


cells appeared morphologically normal, 


When ejected from the micropipette tips, cells returned quickly to their original shapes. It 
noted however, that if larger suction pressures were used (about 5-7cmH,O for a 2um diameter 
pipette) cells were crenated or deformed after release. Ne vertheless, the measured cell volun 


was not dependent on suction pressure for freshly suspended cells if this was n aised above 


2cmH.O, Nor was volume significantly altered if cells were attached to the pipette for periods 
) 


to one minute. Values obtained for cell surface area were not affected by suction pressur 


of attachment at any of the levels mentioned above. 
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SURFACE AREA & VOLUME 


Volume, surface area and diameter were all smaller for older cells and these differences were 
more marked for the rabbit blood samples. However, the cell surface area/volume ratio did not 
decrease with age thus precluding the possibility that erythrocytes become sphered during this 
process. Moreover spherical cells of the diameter measured in this study would have much greater 
volumes than those found here. Nevertheless it was noted that a relatively large proportion of the 
oldest cells were cup shaped stomatocytes with deep depressions inthem. These cells can appear 
roughly spherical but in fact have large surface area/volume ratios. 


It has been shown previously that membrane lipid is lost whilst red cells age (1, . The 
data presented here also suggest some loss of membrane material. However, cell fragmentation 
(10, 15) cannot account totally for the changes in surface area and volume since in that case surface 
area/volume ratio would also decrease (9). It has also been shown that membrane protein is not 
lost during ageing (12) and that cell haemoglobin content remains constant (2, 12). Thus any 
fragments lost during ageing would essentially have to be lipid vesicles devoid of haemoglobin. 
Protein contraction (24) might then account for the decrease in membrane area. Further loss of 
cell water and electrolytes, causing decreases in cell volume without membrane loss, would also 


need to occur to explain the values for the surface area/volume ratio observed in this study. 


With respect to erythrocyte deformability, the reported decreases with age (10, 11) cannot be 
linked to geometric factors. The cell surface area/volume ratio does not decrease with age and for 
the values obtained in this study, the minimum cylindrical diameter which the erythrocytes could 
achieve is slightly smaller for the older celis. That is, for example, the minimum diameter for a 
cylindrical tube through which the upper sample (youngest) of rabbit cells could pass is 3, 05um 
whereas the oldest cells could on average pass through a tube of diameter 2. 94m without their 
volume or surface area having to be altered. On the other hand, the decrease in cell volume and 
concurrent increase in haemoglobin concentration lead to a large increase in internal viscosity as 
the cells age (13) and this will slow the rate at which the cells can be deformed, as shown by Shiga 
etal. (11). Another property which could affect erythrocyte deformability and which might alter 
with age is the cell membrane elasticity. It has previously been demonstrated that membrane 
fluidity and lipid/protein ratio decrease with age (12) and we hope shortly to report results of 


relevant studies of membrane elasticity carried out in this laboratory. 


In summary, we have shown that whereas surface area and volume decrease with increasing cell 


age, the ratio of these parameters remains approximately constant. Thus red cells do not become 


spherical as they approach senesence and their decreased deformability cannot be explained by a 


change in shape. 
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